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IHTR0DUCTI0N 

At a recent conference of the American Geophysical Union it was 
concluded that existing data are inadequate to provide def.nitive models 
of the earth's magnetic field ( Cain , 1975). One of the main reasons is 
that complete spatial coverage of only the field magnitude at any instant 
of time is insufficient to uniquely define a vector field. This problem 
was first recognized by G. Backus and as discussed by Stern 1975, leads 
to large errors in certain terms of the spherical harmonic expansion of 
the earth's field. 

As a consequence of the forgoing problem many investigators have sug- 
gested that a satellite survey of the vector magnetic field is required. The 
p»*ogram proposed by NASA to accomplish this survey has been called f^AGSAT. 

At the present tiire, several study contracts have been let by NASA to 
examine the feasibility of the proposed survey. These studies include 
three different topics; A stable boom design study. An attitude transfer 
design study, and A vector magnetometer design study. 

This report is devoted to one aspect of the vector magnetometer 
design study, a procedure for calibration of the vector magnetoineter. This 
report is numbers in a series of three reports under preparation by the 

University of California in Los Angeles on contract #MAS 5-23660. The preceding 
reports McLeod, 1976, 1977, described a triaxial fluxgate magnetometer 
suitable for the prqjosed mission. 

The body of this report is divided into three main sections. The first 
describes necessary preparations of the test facility. The second describes 
the calibration procedure for individual sensors. The third describes 
procedures for calibrating the sensor assembly. 
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SPECIFICATIONS OF THE VECTOR MAGNETOMETER 

The specifications for the vector magnetometer described in RFP-566233/326 
are included here as Appendix Al, These specifications may be sunmarized 
as follows. The design goal for the magnetometer is to measure three 
cai^onents of the magnetic field in sensor coordinates with a resolution of 
one gamma and an absolute accuracy better than five gamma in a ^64,000<' 
field. This is to be accomplished fpr at least one year over a broad range 
of sensor and electronics temperatures. The magnetometer must be 
constructed so that its accuracy is not impaired by extreme thermal shock 
and exposure to acceleration. The sensor offsets iraist be less than 0.4 y and 
cannot drift nrare than this amount in one year. Noise in the sensor 
outputs must be less than 0.1/ zero to peak in the bandwidth 0.1 to 25 Hertz. 
Magnetic axes of each sensor must be orthogonal within 0.1 degree, known 
to one arc second and stable to five arc seconds. 

COMMENTS ON SPECIFICATIONS 

As pointed out by the author in his proposal for this design study, 
these specifications far exceed those of any magnetometer previously 
carried by satellite and in fact exceed those of most ground magnetic 
observatories. The only detailed report of an observatory which meets 
these specifications is by Yanagihara, et al ., 1973, describing the Kakioka 
observatory in Japan. It should be noted that this observatory has a staff 
of nearly 50; its instruments are on fixed granite pillars in temperature 
controlled rooms; calibrations of all instruments are carried out on an 
almost daily basis. 

The author believes that most of the magnetometer specifications can 
be met by careful electronic and mechanical design. The most serious problem 



hoviever, is angular calibration and stability of the sensors. For example, 
according to Yanagihara et al. , 1973, ordinary magnetic theodolites can 
determine the direction of an unknown field to no better than 6 arc seconds. 

A specially built theodolite in use at Kakioka (A-56, universal standard 
magnetometer) could only obtain 3 arc second accuracy. A more recent 
version of this instrument, the DI-72 has obtained an accuracy of one arc 
second. 

It goes without saying that the accuracy of a calibration can be no 
better than that of the test facility. Thus, from the preceding discussions 
it would appear that the calibration of the orientation of the test coil 
will be limited to about 6" unless special equipment is available. This 
in turn would impose a similar error on the magnetometer calibration. 

In this report we develop an alternative procedure for angular calibration 
that does not requi»-e a magnetic theodolite. Instead, we use measurements 
made by the sensor under calibration to simultaneously determine the 
orier^-ition of both the test coil and sensor. As we will show, this 
procedure is limited by the accuracy with which the magnitude of the 
calibration field can be measured and by the accuracy of the sensor's 
sensitivity and offset. 

DESCRIPTION OF TEST FACILITY 

To meet the specifications of AppendixAl it will be necessary to 
carry out the magnetometer calibration in a well calibrated magnetic 
test facility. At the present time there are only two such facilities 
in the U.S., one at Goddard Space Flight Center and one at the Ames 
Research Laboratory. The Ames facility has a rather small set of calibration 
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coils and no provisions for thermal -vacuum calibration, consequently, 
we will assume that the test will be carried out at GSFC. 

A brief sunanary of some of the facilities available at GSFC for 
magnetometer testing is given in a report by C.A. Harris. 1971. The 
magnetic field component test facility consists of a 22 foot diameter 
three component coil system with remotely isolated magnetometer and 
control instrumentation buildings. The orthogonal field cancelling coils 
are sufficiently large that nearly any magnetic field may be produced in 
a sphere of diameter 3 feet. The main winding of each roil is connected 
to one axis of a 3 axis resonance magnetometer which senses the variations 
in the earth's field and feeds back a current which cancels this field. 

The main winding is also connected to a D.C. field generator which can 
produce fields of up to 60,000J in 0.1;< steps. A second winding on each 
coil is used to cancel the temperature dependence of the main winding. 

A third winding is used to mimimize field gradients over the test 
volume. 

The orthogonal coils are oriented with X horizontal towards magnetic 
north, Y horizontal to the east, and Z vertically downward. The field 
generators can produce 60,000t in Z, 25,000j in X and 6,0002' in the east- 
west direction. Accuracy of field nulling is of order 0.2). The gradients 
across the 3 foot test volume are such that the field does not depart 
from the center value by more than 0.6t. 

DESCRIPTION OF TEST EQUIPMEN T 

Several pieces of test equipment of high precision and accyacy are 
required to carry out the calibrations of the test facility and magnetometer. 
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In this report we assume the following items are available. 

1) 6 digit - digital voltmeter 

2) 3 precision resistors 

3) Proton precession magnetometer 

4} 2-3 component station magnetometers 

5) 2 - (H>tical theodolites (autocollimaters) 

6) 2 - Precision levels 

7) Optical octagon 

8) Brass fixture for rotating sensor assembly 

9) Programmable digital data logging system with analog and 
digital inputs as veil as keytward entry 

Calibration of Test Facility 

To accomplish the absolute accuracy required by the mission requirements 
it is necessary to calibrate the test facility first. Steps in this 
procedure are described in the following subsections. 

^teasure Calibration Coil Constants 

Place a precision resistor in the current loop from the DC field 
generator to the main winding of each coil- Monitor the voltage drop 
across the resistor with the digital voltmeter. Fix the proton precession 
magnetometer at the center of the test volume. Increment the field 
generator in 5000 steps across the full range of the proton precession 
magnetometer 20,000-60,0003). For each step record the field measured 
by the proton magnetometer 10 times and average. Also record the reading 
of the digital voltmeter and convert to current using the known resistance. 
Repeat for the same range of negative field values. 
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Fit a straight line to the field versus current data determining 
the coil constant (slope) and coil offset (intercept). B = kl + B^. 

Also determine the probable errors in these constants using the known 
accuracy of the precision resistor and digital voltmeter as well as the 
standard deviation of the proton magnetometer n^asurements. Assign a 
probable error to any test field calculated on a basis of this formula. 

If the coil offset is not zero or there are systematic departures of 
the field from the linear relationship then a more elaborate calibration 
will be required. This possibility exists since the calibration field 
generator output is nixed with the output of the resonance magnetoir^ter 
i^ich nulls the earth's field. 

Minimize CalibrationrCoil Gradients 

Set the maximum possible field in a coil. Use the proton precession 
magnetometer to n«asure the field along the coil axis at 6 inch intervals. 
Fix the magnetometer at the point of maximum deviation from the center 
value. Alter the current in the gradient adjustment coils to reduce this 
deviation as much as possible. Repeat the survey along the coil axis. 
Again, fix the magnetometer at point of maximum deviation and again adjust. 
Iterate this procedure until the field is as unifonii as possible. 

Test Earth Field Nulling System 

The earth's field nulling system is a servo loop with the sensors 
physically separated from the region in which the field is being nulled. 
Inevitably there are small differences in field bet^en these two locations 
whicii the sensors cannot measure. Because of this there will be a small. 
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residual field at the center of the coil system. This field will have both 
constant and time varying c(»nponents. The magnitude of these con^onents 
limits the accuracy of any magnetometer calibrations. The residual field 
can be measured in the following manner. 

Place a calibrated three conponent fluxgate magnetometer at the center 
of the coil system. Align the sensors roughly along the coll axis. 

Record the output of this magnetometer with the digital data acquisition 
system. 

Perform three series of measurements. First, record about 1 minute 
of data at a sample rate exceeding 120 samples per second. Next record one 
hour of data at a san^le rate of 2 samples per second. Finally record 8 
hours of data at a sample rate of one sample per 5 seconds. Repeat this 
series of measurenffints for the following conditions: Midweek workday, 
midweek night shift, weekend midday, weekend night shift. 

The data gathered in the foregoing experiments should then be subjected 
to power spectral analysis. For each set of conditions separate auto 
spectra are calculated for data at each sample rate. These auto spectra 
can then be plotted as log power versus log frequency on the same 
graph. The data define the magnetic noise in the test facility from 
about 0.3 mini hertz to 100 Hertz. A separate graph for each set of 
conditions show how this noise is a function of time during the work week. 

To define the expected error in any particular magnetic field measure- 
ment we integrate under the appropriate noise spectrum. For a lower frequency 
limit use the reciprocal of the duration of the measurement (the time during 
which we must assume the field is constant). For an upper limit use the 
Nyquist frequency of the magnetometer making the field measurement. The 
rms field error is the square root of the area under the noise power curve 
between the specified limits. 
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Definition of Geographic Coordinates 

Test procedures described below require the definition of an accurate 
geographic coordinate system. This can be done as described below. 

A test table is Installed in the center of the coil system with its 
surface about six inches below the center line of the horizontal coils. 

The table should be made of either granite, marble or glass. It must be 
nonmagnetic and rigidly attached to pylons sunk to bedrock. It should 
have provision for adjusting the surface using precision levels to be 
exactly level. The table should also provide a means for rigidly 
attaching a straight edge to the surface of the table. (This is used 
where repetitive measurements with 180“ rotation are required.) 

Once the test table is installed and leveled the geographic coordinate 
system is established with the use of two theodolites as shown in figure 1. 
The first theodolite is placed north of the test table on a strut or pylon 
sunk to bedrock. This should not be one of the fixtures holding the 
calibration coils. Later it idll be important to determine the orientation 
of the coil axis in geographic coordinates. This can not be done if the 
theodolites are attached to the coil supports. 

The theodolite is leveled by placing a mirror on the test table at 
the exact center of the test coils. The mirror is gradually rotated (with 
tangent screws) until the reflected beam is observed in the theodolite. 

The theodolite is translated vertically and horizontally until it is 
level with the telescope in the same plane as the mirror. 

A second theodolite is placed at the same distance as the first but to 
the east of the coil center. The mirror is replaced by a right angle 
reflector. The second theodolite is translated until a beam from the first 
enters the second theodolite when it is exactly level. 
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The geographic coordinate system defined in this manner has the X 
and Y axes in the hori 2 ontal plane. The X axis is roughly aligned with 
the north coil axis and the Y axis with the east coil axis. The Z axis 
is vertically downward and can be defined by levels when necessary. 

The magnetic field in geographic coordinates can be written as a 
vector sum 

B(GEO) ■ rJjj + By Hy + B^ (1 ) 

In this formula B^. is the magnetic field produced by the i^^ test coil 
and given by the relation 

B. = k. I. 

where is the current and k. the coil constant for the i^^ coil. Also, 

At- 

is a unit vector parallel to the magnetic axis of the i^^coil. In 
geographic coordinates we may write n. as a column vector 

^i " t^Xi X + Hy. Y + Z) (2) 

A A A 

where X, Y, 2 are the orthogonal unit vectors of geographic coordinates 
defined by the theodolites and the vertical direction. Substituting for 
n^. in equation (1) v/e find 

A 

B(GEO) = + Hj^y By + ^ 

A 

+ (riyjj + Oyy By + H y^ B^ ) Y 

+ (n^j^ By, + n^Y By + B^) Z 

This may be written as 

B(GEO) = (n) B(COIL) (3) 

Here B(COIL) is a vector with components made up of the fields generated 

A 

by each coil. The matrix (n) has columns made up of the unit vectors 

A A 

Hy, n^* 
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Definition of Sensor Coordinates 

A second coordinate system fixed with respect to the magnetometer 
sensors must also be defined. In our procedures this system is defined 
by an optical polygon rigidly attached to the orthogonal array of sensors. 
For reasons discussed later this should be an octagon with all sides and 
top face silvered. The octagon Is glued to the sensor assembly with ore 
face as nearly perpendicular to the X sensor axis as possible. Also the 
top of the polygon Is aligned with Its face as nearly normal to the axis 
of the Z sensor as possible. 

For purposes of discussions we assume the sensor assembly consists 
of a cube approximately 3 Inches on a side (typical ring core fluxgate 
array). In this case, the polygon might be mounted on top rf the assembly 
and appear as shown schematicclly In Figure 2 . 

With Sensor Coordinates defined as above the X sensor is approximately 
aligned with the X axis of sensor coordinates and similarly the Y and 
Z sensors are nearly aligned with Y and Z coordinate axes. 

In sensor coordinates we can write the magnetic field measured by the 
.th 

1 sensor as 


Mj. = B-mj j = 1, 2, 3 (4) 

A 

v/here B = B(SEN is the actual magnetic field and m. is unit vector along the 

^ J 

sensor axis, both in Cartesian sensor coordinates. If we write m. as a 

. J 

column vector we have 


A 



"xj "yj ^ * "zj " 


(5) 


where m. . 

* J 

to the i^^ 


A A A 

x-*m. is the direction cosine of the unit vector m. with respect 
1 J j 

axis of sensor coordinates. 


B = B., X + By y + B^ z 


Similarly, we write 
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Substituting in equations (4) we obtain 



= m 

B 

+ 


B„ 

+ 

®, 


• XX 

X 


yx 

y 


2X 

N 

« m 

B 

+ 


B„ 

+ 

m 

y 

xy 

X 


yy 

y 


zy 

N 

= m , 






m 

2 

X2 

X 


yz 

y 
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This nay be written in matrix form as 

N = ( J) B(SEN) (6) 


Tlw matrix (ii) is a matric whose colums are the unit vectors along the 
magnetic axes of the three sensors. The vector H is a vector with con^onents 
given by the actual measurements made by the three sensors. The inverse of 
eq. (6) allows one to obtain a Cartesian vector from the actual measuren^ents , i.e. 

B(SEN) = ( J)'’ H (6A) 


It should be noted that the vector M is not the representation of 
the actual magnetic field in a non-Cartesian coordinate system aligned with 
tte sensor axes. If we want this representation we must write 


B - B m„ B ' M B m 
D ^ y y 2 2 


* “x' [”xx » * ^x > * "xx 


®y‘ C“xy * * "yy * "w 


=B'Cro x+m y+m 
2 *■ xz yz ^ ; 


22 


Rearranging terms this beconies 


4k 

Z] 

A 

2 ] 

A 

2] 


B(SEK) = M B(HEA) 



» t- 
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t^ere B(lffA) * (B„*. B *, B *) is the representation of 6 in the non- 

^ A Jr A ^ 

Cartesian coordinate systen ali^fied with the magnetic axes of the s^sors. 
Using the result, eq. (6) we find 

n = (m^ii) B(MEA) 

or finally 

B(MEA) = H 


Preparation of Ma<yieton«ter Test Fixture 

The test procedures discussed below require ttiat the sensor ass«rt)1y 
te oriented in a variety, of precisely known orientations relative to 
geographic coordinates. This is accoiqilistod by rigidly attaching an 
optical polygon to the sensor ass^bly. Reflections frm the faces of 
ttose polygons are aionitored by the two theodolites and used to calculate 
the precise orientation of the polygon. The difficulty with this procedure 
is the limited field of view of the theodolites. It is not possible to 
perform arbitrary reorientations of the sensor asseiri>ly and have a face 
of the polygon nearly normal (within about a half degree) of the theodolite 
optical axis. 

To solve the foregoing problem we utilize a magnetiXReter test fixture 
which nukes possible reasonably accurate rotations about two orthogonal 
axes. An example of such a fixture is shown in Figure 3. The device 
shown is an earth inductor or magnetic theodolite. Its purpose is the 
accurate determination of the direction of an unknown magnetic field. 

It utilizes a spinning search coil to indicate vdien the search coil 


rotaticm axis is aligned with the ambient field. The azimuth and elevatiwi 
of the rotation axis are read from horizontal and vertical circles. 

We propose to modify such a (tevice as shmn schematically in figure 4 
The search coil drive assembly is replaced by holloM tubes ccmcentric 
with the elevation axel (horizontal axis). This allows unobstructed 
observaticm of one face of the (^tical octagon when moimted within the 
fixture. The search coil is replaced by a mounting plate with attacNmnt 
cl«93. Hie platform and clamps are provided with adjustmmts idiich liable 
the experimentor to align the optical axes of the polygon with the rotaticm 
axes of the fixture. 

The procedure for setting up such a fixture would te as described 
telow. Place the magnetometer sensor assembly on the mounting plate of 
the fixture. The plate shcxild be designed such that the center of the 
optical polygon is close to the intersection of the vertical and trarizontal 
axes of the fixture. Place the fixture <mi the test table with the 
horizontal axis pointing north and the center of cme face of the octagon 
in line with the north theodolite. Move the fixture along the north- 
south line until the center of an orthogonal face is roughly aligned wite 
the east theodolite. Level the fixture using the three leveling screws 
and levels attached to the base of the fixture. Rotate the fixture in 
azimuth with a tangent screw until the normal to the northward face lies 
in a vertical plane passing through the north theodolite. Next, use the 
adjustment screws on the mounting plate to bring the normal to the east 
face into a vertical plane passing through the east theodolite. Repeat 
these steps until both fa:es of the polygon are nearly orthogonal to the 
two theodolites. When thi.' is achieved the top surface of the polygon should 
be almost level. 
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Nom continue with fim adjust^nt of fixture level to bring the 
azinuth axel (vertical axis) vertical. Do this by rotating in aziaiuth 
by 45** increcaents. tf the successive faces do not produce centered ioages 
In the theodolites thst azimuth rotation axis is not truly vertical. Hake 
it so by fine adjustments of fi»*-ure level. 

When the foregoing procedure is conpleted it should be possible to 
perform rotations through angles of integral multiples of 45** and maintain 
reflected images of crosshairs within the field of view of the theodolites. 
Because of the inaccuracies associated with the bearings of the test 
fixture it Is not expected tt»t the rotaticm axis will be stable to much 
better than same fractitwi of a minute of arc. Also errors associated with 
the verniers on the azimuth ami elevation circles will limit the accuracy 
of rotatitxi angles to some fracticm of a minute. 

This problem is not important. Actual alignasnt of the sensor coordi- 
nates in geographic coordinates is determined by the theodolite reflections 
from tte precisely constructed faces of ^e optical octagmi. The primary 
purpose of the fixture is to provide sufficiently accurate rotations that 
the reflected images renain within the field of view of the theodolites. 

Relation between Sensor and teographic Coordinates 

In the previous section we described a procedure for setting up a 
test fixture which holds the array of sensors to be calibrated. The 
sensor assembly has attached to it an optical octagon with its top and 
side faces silvered. The faces of this octagon must be constructed so 
that they meet at 45** angle and are flat to better than one arc second. 

The top must be at 90** to all faces within the same accuracy. The fixture 
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IS sufficiently precise that it can rotate the sensor assembly through 
angles in increments of 45® so that images reflected to the theodolites 
viewing the orthogonal faces of the octagon I'emain in the field of view. 

It is not expected that these images will be exactly centered. The exact 
orientation is thus determined by using the theodolites as eutocol lima tors 
reading the precise orientation to a fraction of an arc second. 

As an exan^le, consider the case shown in the top of Figure 5 for which 
the X and Y axes of sensor and geographic coordinates are approximately 
aligned. In actual fact, the X axes are not exactly coincident as shown 
in the lower left. The azimuthal angle ^ between the X axes in the 
horizontal plane is measured by the horizontal displacensnt of the image 
in the north theodolite. Similarly the polar angle is given by the 
vertical displacement of the image in this theodolite. Using these angles 
wa can express the direction cosines of the X axis of sensor coordinates 
in geographic coordinates as 

X^(GEO) = [cos^jjSinrjj, Sin^^Sin^^^, Cos^ (7) 

A 

Similar measurements made with the east theodolite give the Y axis as 

YjCGEO) = [Cos{^y+ 90) Sin^y, Sin(^Y+90) Sin9Y.Cos9Y) (8) 

A 

Finally the direction cosines of the Z axis are found by the requirement 

A A A 

that Zj is orthogonal to X^ and Y^. 

Z^(GEO) = X 3 X Y 3 (9) 

/V /N 

Where X^ and Y^ are given as above. 

The transformation from geographic to sensor coordinates is found by 
using the rule that the rows of the transformation matrix are the unit 
vectors of the new coordinate system expressed in the old system. Thus 
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B(SEN) = (R) B{GEO) (10) 



Successive rows of (R) are given by the el«ients of the vectors appearing 
in equations 7, 8, and 9 respectively. 

With the results presented above and in previous sections we can 
easily relate the measured magnetic field to the field produced by che test 
substituting equation (3) into (10) we have 

B(SEN) = (R)(n)(BC0IL) 

then, substituting this result in equation (6) we find 

M = (J)(R)(n)B(C0lL) (12) 

Equation (12) graphically illustrates the fact that any measurement 
made in the calibration facility couples the unknown direction cosines 
of both the coil and the magnetometer axes. Also it is apparent that this 
equation is nonlinear in the unknowns. If we perform a number of experiments 
creating a set of such equations we must solve a set of simultaneous, non- 
linear equations. This problem is further complicated by the fact that 
each column of (y) and (iq) is a unit vector, i.e. there are only 2 unknowns 
rather than three per column. Thus, altogether there are 12 unknowns to 
be determined experimentally. A procedure for accomplishing this is 
described in a subsequent section. 
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CALIBRATION OF INDIVIDUAL SENSORS 

The second major step of the vector magnetometer calibration procedure 
is the calibration of individual sensors. This calibration includes the 
determination of sensor sensitivity offset, noise and drift. In addition, 
it involves the measurement of the effects of temperature and magnetic 
fields orthogonal to the sensor. In following subsections we discuss each 
of these procedures in detail. 

Measurement of Sensor Sensitivity and Nonlinearity 

The component of magnetic field parallel to the axis of a linear 
fluxgate sensor may be written 

®i " *'i^i ^ °i ^^3) 

In this equation is the sensor sensitivity in gairana/volt, and 0- 
is the sensor offset in gamma. To determine these constants we must 
apply known magnetic fields B and measure the sensor output voltage, 

Plotting the applied field as a function of output voltage, the slope of 
a best fit straight line determines the sensitivity and the intercept 
determines the offset, 0^. Systematic deviations of the measurements from 
the best fit line indicate that the sensor is not truly linear. 

Normally, the sensitivity of a magnetometer is determined from 
only tv/o sets of measurements, one at zero field and one near full range. 

In the following paragraphs we show this procedure is not sufficiently 
accurate to meet the requirement of 1 gamma absolute accuracy over the 
full dynamic range of the sensor. Consequently, it is necessary to make 
a large number of measurements across the full dynamic range of the sensor. 
These measurements are then fit by a least square straight line. This 



19 


procedure Improves the accuracy of the sensitivity determination and has 
the additional advantage that a plot may be made of the deviation of the 
observations from the best fit. 

To estimate the accuracy of a particular magnetic field measurement 
based on equation (13) we note 

= S(k V) +^(0) 
but 

&(k V) ^ rli + iV-, 
kV k V-* 

iB = kV + £X] . /(o) (,6) 


A crude estimate of the values of k and 0 and their errors can be made by 
using only tv/o sets of measurements. Suppose Bq and Vq correspond to zero 
field and B^ and correspond to 50,0003 and 10 volts. Then, from 
equation (13) 


k = 



and 0 





( 17 ) 


We can show 


and 


£JS. - 2 {1 . 2 |V 

k - B, V, 


T " Csf " ‘ <V7> “ 


Since B, 0 we have finally 


(18) 


(19) 


iB r 

B 


kV[^ + 


3 51 -, 

Vl 




] B, 


For a voltage V close to V.j, kV so that 


5_B ^ 3 5 B , 
■ B By 


4 5V 

h 


(20) 


( 21 ) 
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For most test facilities it is difficult if not inpossible, to generate a 
calibration field accurate to better than IK in 50,00GJ'. Also available 
digital voltmeters can measure to an absolute accuracy of about 50/f V in 
10 volts. Thus we take 

^B/B = 11750.000K' = 20 X lO"® 

SV/V = 50 X 10“^/10 = 5 X 10~® 

It is evident that ♦ accuracy of the calibration field is most important 
in determining the final accuracy of the calibrated sensor. Substituting 
in equation (21) we find for full range 

“ 80 X '0-« or <B^^=r 

The minimum error occurs for zero field and depends on the error in the 
offset 

(^B)^i^ = S0= [WBj + SV/V, 3 B, 

Nun»rically 


These errors are unacceptably large and demonstrate the need for a 
more accurate calibration procedure. In the following subsection we 
describe an independent procedure for determining offset to an accuracy 
of ^^O.lX. Next, however, we consider the inprovement in accuracy obtained 
by the least square procedure mentioned above. 

In fitting a straight line to the calibration date we write 

®v " \ ^ V = 0, 1, 2,...N 

Here V is in index referencing one of the pairs of N+1 calibration measurements. 
In this procedure we choose as the dependent variable, because it can 
be measured more accurately than In the usual least square error 
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analysis it is assumed that the dependent variable is error free. Using 
results derived in Hildebrand, 1956, we can show the coefficients Aj^ 
satisfy the set of normal equations 


(c) A V 


where 



(23) 


The error in each coefficient is given by 
2 


where 


jfl = fc'^l f^^l 5^ 

\ ^kk ^N-1 ^ " RMS 




(24) 


Note ^j^vjs square deviation of the observations from the predic- 

tions; i.e. fS in our previous discussion. 

The normal matrix c is easily inverted so that 


-fv, (N*l), 


and 


N 5 N 9 
|c| = (N+1) ^ V/ - ( £ V.)'^ 

i=0 ^ i=0 ^ 


To estimate the errors we must assume some model of the calibration procedure. 
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Therefore, assume 10 volts corresponds to 50,000) and that we step through 
the dynamic range 0 - 50,000j in 500^ steps, i.e. 100 sters of 0,1 volt. 
Then (N+1) = 101 


N 100 , 

£ Y. = .1 ^ (i) = .l[y(100)(101)] = 505 

1=0 ^ i=0 ^ 

N p 9 100 „ 

£ V. = (.1) ^ (i)"^ = ,01[338,350] = 3333.50 

i=0 ’ i=0 


Thus |c| = (101) (3383. 50) - (505)^ = 88,708.5 

^.v.2 

Hence, (c’^)qq = = *03902 

(c‘^)ll - ^ = .0011648 
For N = 100 v/e have 


Thus 


SIKqX .2 


Sms 


(25) 


This result should be compared to our previous error estimates based 
on two pairs of measurements. Using equations (18) and (19) and ignoring 
SV we found 




The use of 100 pairs of measurements should improve the calibration by about 
a factor of 5. 



23 


The error in a field measurement was according to equation (16), 

= kY [^ + ^] + ^0 
since JAq = So and we have 

‘ ^ * 5 * 'O'®] " -2^15 

^RMS " 

While the foregoing error is quite acceptable it depends on the 
linearity of the magnetometer. If the observations systematically depart 
from a straight line ^dmc will be larger than 13T, and SB will be 

KrlJ nioX 

proportionally larger. In this cace, it might be necessary to use higher 
order functions to fit the observations. 

On a basis of the preceding analysis we recommend the procedure 
shown schematically in figure 6, to determine sensitivity and offset of 
each sensor. Attach the sensor to a test fixture which has provisions for 
slight rotations of the sensor about two axes. Place the fixture on the 
test tabie with one edge against a north-south straight edge attached to 
the table. Apply maximum field (-50,000j) in the north direction. 

Rotate the sensor around a vertical axis to obtain maximum output voltage. 
Next, rotate the sensor about a horizontal axis again maximizing the 
output voltage. Repeat these steps several times until the best possible 
alignment of sensor magnetic axis and calibration field is obtained. 

Once alignment is achieved perform a sequence of measurements 
determining sensor output voltage as a function of calibration magnetic 
field. To determine the precise input field use a digital voltmeter to 
monitor the voltage drop across a precision resistor placed in the field 
generation drive circuit. Using the calibration coil constant determined 
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by th“ method described above and the measured coil current 
(Icon "" V/R) calculate the applied field. Also use the digital voltmeter 
to measure the sensor output voltage corresponding to the input field. 

To obtain sufficient accuracy perform 101 pairs of measurements with 
the input field incremented in 1000 steps over the range -50,000 to 
+50.000J'. Enter the table of measurements, versus into a computer 
program which fits a least square line to the data determining k, 0 and 
their associated errors. The program sho. d also calculate and plot the 
deviation of the predicted field from the observed field. If the resulting 
time series is Gaussian with zero mean the sensor is linear. 

This procedure is repeated twice more for the remaining two sensors. 
Together these three experiments define the sensitivities and offsets 
required to measure the three components of the magnetic field provided 
the orientation of each sensor is known in inertial space. 
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Measurement of Sensor Offset in Low Field 

An independent, and more sensitive determination jf sensor offset 
can be made by tests performed in a low field environment. This procedure 
is somewhat easier than the one described in the preceding subsection 
and because it does not require a large test facility, we recommend 
that it be used for long term monitoring of temporal drifts in offset. 

The sensor offset is defined by equation (13) 

B = kV + 0 

Solving for the magnetometer output voltage we find 
'' k 

Suppose we place the sensor in alignment with a weak field Bq, the sensor 
output will be 



Now reverse the direction of the field either by a 180 rotation of the 
sensor or by changing the sense of the current producing the field. The 
output voltage is 



Adding the two measurements and solving for 0 we have 
Y +Y 

0 - (4l~) (26) 

Note the sensitivity must be independently determined to calculate offset. 
The error in 0 is roughly 

io “ t 0 

from equation (25) ^k/k t 7 x 10'°. Hovjover SV/V is of order 1. The 
difference is that Sv should not be the precision of the measurement device, 
but the fluctuation induced by variations in the ambient field and by 



26 


instn^nt noise. In an unshielded, industrial environrent fluctuations in 
voltage JV correspond to field fluctuating of order IJ. Since ring core 
offsets are of tte sanse magnitude we expect SV - V. Clearly accuracy can 
be obtained only by repeated measurements in a shielded environment. 

To carry out a determination of offset we recommend the following. 

Place the sensor on a fixture which allows an approximate 180** rotation. 

(An accuracy of a few degrees in this rotation is sufficient.) Place the 
sensor and fixture inside a set of ctwicentric ou metal cans. Cover 
each can with its i^tal cap -*xcluding the earth's field fran the 
interior of the innermost can. A remanent magnetic field of a few gamma 
magnitude and unkn<^ directim will rerain in the can. Since the direction 
of this field cannot be changed the orientations of the sensor ixist be 
reversed. Perform a series of me'surements of tte magnet(xneter output 
voltage, rotating the sensor 180** before each measurement. Proceeoing 
through the table of measurenents, average pairs of readings and calculate 
offset. Average the offsets so determined to obtain a final, more accurate 
value. 

Measurei?»nt of Sensor Noise in Low Field 

In preceding sections we have discussed offset as if it were a constant 
property of the sensor. In fact, offset usually changes with time in a 
random manner. For fluxgate sensors, a frequency spectrum of these changes 
is typically inversely proportional to frequency. Thus, on a short tine 
scale, variations in offset are quite small. In general, these variations 
are divided into tv<o categories depending on their time scale. Variations 
which take longer than some reference time are called offset while those 
that take less time are called noise. A typical rr nee time is 
between a day and a v/eek. 
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Sensor noise is best charactarized by a power spectrum of the sensor 
output when located in a zero field environment. In such a situation it can 
be assumed that all variations in output are a result of the instrument 
rather than the asdiient field. 

We suppose that the sensor has been placed inside a set of concentric 
mu metal cans as described in the previtMis subsection. The sensor 
output is recorded by a digital data acquisition syst^. Usually, it is 
necessary to amplify the magnetocster output voltage prior to digitizaticm 
or quantizatitm noise created by digitization will dominate the spectrum 
at higher frequencies. Also it is ieportant to low pass filter the sensor 
output volta^ with a time constant twice the sampling interval. This 
eliminates the problem of aliasing noise power from high to low frequencies 
in the digitization procedure. 

Between 1(K)0 and 50G0 sonples of the sensor output should be taken 
and then read into a congniter program for spectral analysis. This program 
estimates the noise power in a frequency band corresponding to about 10/T 
to l/2it idiere T is the duratiw of the series of measuren»nts and At the 
soling interval. For 5000 samples we have a ratio of upper to lower 
frequency limits of 250 or 2.4 decades of frequency. 

To cover a wider band of frequencies in an efficient manner, we must 
repeat the above experiment with progressively higher sarnpling rates. A 
convenient set of experiments is stmciarized schematically in Figure 7. 

The frequency band from 10~^ Hertz (1 day) to 50 Hertz (100 sanples per 
second) can be covered in three experiments. These are one day of recording 
with one minute san^les, one hour of recording with one second san^les and 
one minute of recording with 100 san^les per second. 

In order to consider the effects of quantization noise in this 
measuren^nt we must assume son» typical noise spectrum. For reference in 
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Figure 7 ve have plotted a spectruo obtained at UCLA of a ring core 
fluxgate sensor scaled to ” 8000V. Ffeasurecents made with ring core sensors 
scaled to ■ 64,(KX)J are not significantly larger. For the curve sh<^ the 

_3 

nns noise power over any three decades of frequency is 26 x 10 

Quantization noise has a flat (idiite) spectrum with magnitude 


Pq{f) • 


12 fg 


where is the quantization level of the input data and f^ in the bandwidth 
of the measurraeots. Using the Byquist fr^uency as the bandwidth we obtain 

Pq(f) = (27) 

According to this formula must satisfy 



If we wish the quantization noise to be below the expected instrunmnt noise 
in each experim^t (see dashed lines in Figure 7} then va raist have 
SB .V For an instnm^nt scaled to + b0,000 this corresponds to 

SB/B = 2 X 10'^ 

22 

or 1 part in 5 million. Equivalently this is about 1/2 . 

ir suppose instead the quantization -corresponds to 1 gasisna, then 
Pp = At/6. This is plotted as a dashed line for the experiment At = .01 sec. 
Quantization noise is then several orders of magnitude greater than expected 
sensor noise. 

From the foregoing discussion it is apparent that the magnetoireter 

22 

signal irust be quantized to about 1 part in 2 if we are to observe 
instrument noise when the dynamic range is + 64,000>. If the sensor has 
only a digital output quantized to 1^' in + 64,000T, sensor noise cannot 
be measured. If the sensor has an analog output, its output signal must 
be amiplified so that the least significant bit of the data acquisition system 
correspond to about .02 V. 
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For exafl|}1e, suppose, 50,(KK)][ corresponds to 10 volts. Then 
.023' = K/t volts. If the data acquisition system provides a 16 bit converter 
over a range 0-10 volts, the least significant bit is ^150/ volts. An 
an^lification of about 40 is required to fneasure the sensor noise. 

Our recommended procedure for cieasuring the noise of the sensors 
is suas»rized as follows. Place the sensor in a shielded container. Low 
pass filter and saaple the (Hitput signal. Record the san4)les on digital 
tape and subject the data to spectral analysis. Plot the results as log 
power versus log frequency. Three separate experiments shtnild be perfonaed. 
These corresp<wid to one day of 69 second sac^les, one hour of one second 
sasples and 1 minute of .01 sairples. To be ^aningful the effective 
quantization of the data should te about .02V in 50,(X)0t or 1 part in 2^^. 

Kgasure?rent of Teesoral Drift in Sensor Offset 

In the preceding section we pointed out that sensor offset is time 
varying. By definition, changes in offset on a time scale longer than a 
day are called ten^oral drift . The drift can only be measured by repeated 
measurements at widely separated tin^s. Since it iinxild be difficult to 
routinely carry out such measurements in a large test facility we recommend 
use of the low field offset determination procedure. 

To determine the nature and magnitude of tcn^oral drift proceed as 
follows. Once a week place a continuously operating three component magneto- 
meter in the test fixture inside a shielded can. Carry out a series of 180® 
rotations and calculate offset. Repeat for the reir.aining two axes. Return 
the magnetometer to an isolated location where it continues to monitor the 
earth's field. Plot the three sensor offsets as a function of time for about 
one year. At the end of the year determine the naan offset and rms deviation 


about the mean. 
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Keasurerent of Teciperature Dapendence of Sensitivity and Offset 

It is Mell established that the sensitivity as %«e11 as offset of flux- 
gate nagnetometers is a function of teH^)erature. Since the sensor electronics 
is separated froa the sensor assesi^ly, both the electronics teirperature 
(Tg) and the sensor teaperature (T^) are significant variables. Taking 
this dependence into account %«e rewrite equation (13) as 

B, ■ Tj) V, * 0,(Tj Tj) (2B) 

Considering the two parameters k and 0 as functions of and we 
aake the crude assiaption 

P(Tj Tj) = P(Te„ t aJTj + WTj (29) 

In this expressiOT P is either of the two parameters k or 0. The subscript 
"o" designates the nominal operating temperature (assun» rc(^ temperature) 
of the sensor and electronics. ilT is then the deviation of either tempera- 
ture from the nominal value. Finally, a and b are teaperature coefficients 
for the parameter. 

It should not be ejq)ected that equation (29) applies over the full 
range of operating temperatures for the roagnetoneter. Experience at UCLA 
has shovni that sensitivity is often a quadratic function of both temperatures. 
However, because of the rather large temperature dependence of typical 
fluxgate sensors the MAGSAT roagnetwneter will certainly have thermal 
enclosures about both the sensor assembly and sensor electronics. Within 
this enclosure teaperature variations should be so small that equation (29) 
will be an adequate approximation. 

To determine the temperature dependence of k and 0 we proceed 
as follows. Place the sensor assembly at the center of a three axis 
calibration coil facility. Align the three axes of the sensors as closely 
as possible with the three coil axes. Place thermal enclosures about both 
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the e1ectr<Kiics and the sensor assembly. Using hot and cold gases bring 
both enclosures to nominal operating ter^eratures. Allow sufficient 
time for both units to come to the temperature of their respective enclo- 
sures (approximately one half hour). Determine the sensor offsets and 
sensitivities for each axis by applying a sequence of calibration fields in 
the appropriate calibration coil. Use least square procedures to calculate 
k. and 0.. feeping electronics tatiperature constant increase sensor 
temperature by 5®. Allow the sensor to come to thermal equilibrium. Again 
carry out a series of measurements determining k. and 0^. Next, ino'ement 
sensor teo^erature by 10®, establish thermal equilibrium, and determine k. 
and 0-. Finally, increesnt another 15® and determine k^ and 0-. (Note 
this sequence defines the parameters for = 0®, 5®, 15® 30®.) Now 
return to nominal sensor ten^erature decre.nanting temperature by the same 
amounts as it was previously incrensnted. Continue to negative 
tenparature deviations using the same schenie = -5®, -10®, -15®). 
Finally, return to nominal sensor ten9)erature. 

Plot the parameters k. and 0- as functions of T^. Fit straight 
lines to the three points at = -5®, 0®, +5®. The slope of these lines 
are the coefficients b in equation (29). The zero intercepts are the values 


"'Teo.^So*- 


The foregoing procedure is now repeated holding sensor teinperature 
constant and varying electronics ten?>erature. Again a plot is used to 
determine the coefficient a in equation (29). 

There is a remote possibility that all three sensors might have 
extrema for k- and 0^ at nearly the same teinperature. If this were the 
case it v/ould be desirable to choose these coimon temperatures as the nominal 
operating ten^eratures and repeat the measurement of a and b about these 
new temperatures. It should be sufficient to use only three temperatures. 
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i.c. , 4T = -5®, 0® and +5®. For these points the coefficients a and b will 
be smaller and the magnetometer sensitivity to temperature will be 
reduced: 

Measurement of Orthogonal Field Effects on Sensitivity and Offset 

In the idealized theory of the fluxgate magnetometer the 
sensor response is unaffected by fields orthogonal to the sensor axis. 

Sojne sensors, however, have experimentally shown changes in properties 
when the orthogonal field is very large (McLeod, personal conanunication, 
1976). Proper calibration of the MAGSAT magnetometer should include a 
demonstration that sensor properties do not change as a function of ortho- 
gonal fields. 

To determine the effect of a strong orthogonal field we repeat 
the determination of sensitivity and offset described earlier. In this case 
however, the measurements are taken with a constant 60,000^ field in two 
directions normal to the axis under calibration. Both sensitivity (k) 
and offset (0) are determined as before. If the values for k and 0 differ 
by more than the experimental error a more elaborate calibration is required. 
If this is the case, we recr.sTjend that k and 0 be repeatedly determined 
for different values of orthogonal field. A possible series of measurements 
would start at -60,000j and proceed to +60,000J in 20,9003 steps. A plot 
of k and 0 as a function of orthogonal field should be made. Mote this 
imist be done for both possible orientations of orthogonal field. 

If the suggested effect exists it will hopefully be small. 

In such an event, it may be possible to make a linear approximation of 
the dependence of k and 0 on orthogonal field magnitude. Data would be 
corrected by using tne zero field constants to calculate the fields 
orthogonal to the sensor. These fields would then be used to determine 
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correct values of k and 0 and then a second calculation of the ambient 
field would be made. 
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CALIBRATION OF SENSOR ASSEMBLY 

In this section we consider the most complex aspect of calibrating 
a vector magnetometer; determining the direction cosines of the magnetic 
axis of each sensor in a geometric coordinate system fixed in the sensor 
array. In section we showed 


B(GEO) = (n) B(COIL) 

(3) 

M = (/) B (SEN) 

(6) 

B(SEN) = (R) 8(GEO) 

(10) 


where "(GEO)" indicates a vector in Cartesian geographic coordinates (fixed 
in earth); "(SEN)" indicates a vector in Cartesian sensor coordinates (fixed 
in sensor assembly); "(COIL)" indicates a vector constructed from the magni- 
tudes of the three fields generated by three, nearly orthogonal calibration 
coils, and M is a vector constructed from the magnitudes of the three 
fields measured by three, nearly orthogonal sensors. The matrices (n), (u)» 
and (R) are transformation matrices constructed from unit vectors. Columns of 
(n) are unit vectors of the calibration coil in geographic coordinates; col- 
umns of (u) are unit vectors of the magnetic axis in sensor coordinates; 
columns of (R) are the unit vectors of the geographic coordinate system in 
sensor coordinates. Only the matrix (R) is orthogonal (i.e. R^ = R~^) 
since the sets of unit vectors M- and are not orthogonal, i.e. 

*”i**^j ^ ^ n^»nj ^ 0 

In the'^e expressions we assume that accurate determinations of the coil 
constants and sensor sensitivity and offset have been already performed 
as described in earlier sections. 

The basic problem is to determine the direction cosines of the 

^ A A 

magnetic axes m^, m^, m^ in sensor coordinates. This must be done by 
generating known fields B(COIL) and measuring the output M. 


From 
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(3), (6) and (10) we have 

M = (m)^ (R) (n) B(COIL) (12) 

The transformation (R) from geographic to sensor coordinates can 
be experimentally measured as described in previously. Also (1 and 
B(COIL) are knov«i experin«n tally. Clearly, we cannot find the elements 
of (u) unless we already know the elements of (n) or unless we determine 
them simultaneously. We discuss these two cases separately in following 
sections. 


Determination of Magnetic Axes Orientation Given Direction Cosines of Coil 
System. 

If (n) is known in eq. (12) as v/ell as B(COIL) and M w® 
have three equations and nine unknowns. If we perform tv/o additional 
experiments using different calibration fields we will have nine equations 
and nine unknowns which enables us to solve for the elements of the 
matrix (p). 


The simplist sequence of calibration fields to use in this procedure 
is one in which successive calibration fields are parallel to the 
coil axes. We thus have the three vector equations 


where 


M. 


B^(COIL) 



(R)(n)] B.(COIL) 
B.(COIL) 




i = 1.2,3 



1=j 


We assume for convenience that the geographic and sensor coordinate 
systems have been constructed such that their axes are nearly aligned 
with those of the calibration coils and magnetic sensors respectively. Then 
both (/') and (l) are close to being identity matrices. If in addition we 
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align the sensor and geographic coordinate systems as closely as possible, 
the matrix (R) is also nearly the identity matrix. Thus the expression 
[(u)^ (R)(n)] is also close to (I) and V'/s expect the measured field B^(HEA) 
to have a large consonant in the axis nearly aligned with the calibration 
coil and small components orthogonal to this direction. 

He can combine the results of the three successive measurements 
into a single matrix equation 

(B^^) = Rn] (B^) (30) 


where the three vectors H. form the columns of (B„) and B. (COIL) the 

*1 n -I 

columns of (B^). The matrix (B^) is diagonal by our choice of calibration 
procedure , \ 



The matrix (B ) is approximately diagonal 

^11 8^2 

~ I ®21 ®22 ® 23 1 


B 


V 


21 

31 


B 


B. 


32 ^^33 


'V 


according to eq. (6a) 


B(SEN) = (a;^)'^ M 

so solving eq. (30) for we obtain 

C^"^) - 1 = R (B^B^'b"’ 


(31) 


The matrix (B)^ is diagonal hence its inverse has elements which 
are the reciprocals of the elements of (B^.). Furthermore, post multiplication 
of a matrix (8^^^) by a diagonal matrix is equivalent to m.ultiplying each 
column of (B^^) by the corresponding diagonal element. Thus 
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A 


(b) = 


m c 


-1 


f ®21^®1 ®22^®2 ®23^^3 j 

\® 31^®1 ® 32^®3 ^32^hy 


(32) 


The matrix (b) is constructed by •'ormalizing the three sensor 
measurements in each experiment by the corresponding calibration field. 
Thus, 

( = Rn b*’ (33) 


Errors in Magnetic Axis Orientation Given Direction Cosines of the 
Calibration Coils 

The preceding result suggests that the direction cosines of the 
coils are quite easily measured if the direction cosines of the calibration 
coils are knovvn. The magnetometer array is aligned with the calibration 
coils, three calibration fields are applied in three successive axes, 
and the resulting magnetometer measurements are used to calculate the matrix 
(p) according to eq. (32). A consideration of errors in this procedure, 
however, indicates that this simple procedure must be modified somewhat. 

To examine the errors in (33) we v/rite each of the measured matrices 
as the sum of the true matrix and cn error matrix. 

(R) = (R”) + (/) 

(0 = (n°) + (e'') 

(f') « (b‘’)° + (c") 

thus 

(/)-^ = (R° + t^) (n“ +e"‘) ((b‘^)» + 
or 

= R'’n°(b"b” + [R”e'’(b‘^)° + E^’n°(b"b° + R°n” e^] 
where we neglect terms of second and third order in the error matrices. 


1 


I 


1 
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- R‘>n“{B‘'')" + [c^ + e"' + (34) 

The errors in any element of (y) is approximately the sum of the errors 
in the elements of the measured matrices. Consider first the effect of 
magnetic field measurement errors (e ). The magnetometer output in 
a noisey test site will fluctuate by about + ly. If our calibration 
field magnitude is 50,000y, the error matrix (e ) will have elements 
of order + 1/50,000 = + .00002. Let us suppose this error occurs in 
an element that should be exactly 0.0, i.e. the magnetic axis is exactly 
orthogonal to one of the geometric axes of the sensor assembly. Then 
the angular error in the orientation of the magnetometer axis is 

«e = [Arc Cos (-.00002) - Arc Cos (+.00002)] 
or 

60 = .00114° = 4' 

The calibration error in (y) due to this source alone exceeds 
the design goal of one arc second. We note that if the site and 
instrument are both quiet, quantization error will probably be of this 
order. Note that ly resolution in + 65 Ky requires a 17 bit converter 
on the magnetometer output. 

Errors due to the transformation from geographic to sensor 
coordinates are expected to be small. Large dynamic range autocollinators 
have accuracies of order one arc second (Schneider and Kolany, 1967). 

The largest error is likely to be due to measurement of the 
direction cosines of the calibration coils. For example, a magnetic 
theodolite such as that shown in Figure 3 has an accuracy of about 6 arc 
seconds. More elaborate versions of this instrument described by 
Yanagihara, 1973, have accuracies of 3" and 1". 
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If we suppose we use the best magnetic theodolite to determine 
(n), a wide dynamic range autocollimator to determine (R) and improve 
our magr^'tic field measurement to correspond to one arc second, our 
errors in (u) are still or order 3'*. In a 50,000y field this corresponds 
to about 0.7y. 

To achieve the equivalent of one arc second accuracy in the 

magnetic field measurements of this calibration procedure we require 

6B = 0.2y. But 0.2/2(50,000) = 2 x 10‘® or 1 part in 500,000. This 

19 

is about one bit in 2 . The best available digital voltmeters have an 

absolute accuracy of about 50tiV in 10 volts or 1 part in 200,000. Thus, 
neither a digital magnetometer using a 16 bit converter nor an analog 
magnetometer using the best available digital voltmeter will provide 
sufficient measurement accuracy to give the equivalent of one arc second 
accuracy in the determination of (r). 

From the preceding argument it appears necessary to make 
repeated measurements of the field to obtain the necessary accuracy 
through averaging. However, because of the quantization problem discussed 
above, this may not work. Since instrument noise, and possibly test site 
noise, may both be smaller than the quantization level, the actual 
magnetometer output may nut be uniformly distributed within the quantiza- 
tion interval. In such situations repeated measurements do not increase 
accuracy. 

To eliminate the foregoing problem, and also to improve the 
signal to site noise ratio, we propose the following. All three 
calibration coils are simultaneously driven in phase by a sinusoidal 
current nf precisely known frequency (=0.1 Hz). A large amplitude signal 
(50,0C0y) is used in one axis as before, and small amplitudes are used 
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in the remaining two axes (IOy). The drive currents and magnetoiseter 
outputs are oonitored by a digital data acquisition system for 100 
cycles (20 minutes). The a^litude and phase (should be zero) of all 
signals is determined by least square fitting of a sine wave of the 
known drive frequency to the measured data. 

The small fields used in the two coils orthogonal to the rain 
drive coil will force the magnetometer signals in the corresponding 
axes to cross several quantization levels. In this manner, the small 
fields produced in these axes by the main calibration field can be 
accurately measured because of a uniform distribution of the sensor outputs 
across a quantization level. 

It should be noted that using this procedure the calibration 
uatrix (B ) constructed from the measured anplitudes is no longer 
diagonal. As a consequence, we mast use matrix inethods to determine the 
matrix, 

b = 

With this Kodificaticn the procedure outlined at the beginning of the 
section should be adequate. 
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PeterBiination of Magnetic Axes and Coil Axes SiBailtaneousl y 


If ve assume that the direction cosines of the calibration coils are 
unknown vie can still use eq. (12) to determine the orientation of the magnetic 
axes as well as those of the coils. In this case eq. (12) constitutes three 
equations for IS unknowns. If as before we apply three calibration fields 
in three orthogonal directions we have (eq) 30 which is a set of 9 equations 
for 18 unknowns. If we include the unit vector constraint on each column of 
(^) and (q) we still have 12 unknowns. Clearly there is insufficient 
information to solve for the unknowns. 

Additional information can be obtained by reorienting the sensor assembly 
with respect to the calibration coils. This changes the elefe»nts of (R) 
relating the unknowns to the n^asurements. Using the superscript^ to 
designate the particular orientation of sensor and geographic coordinates, 
eq. (30) can be written 

RrJ = (BjJ) (B*)'^ *<=1,2.... (35) 


For example, using two orientations we should have a set of 18 ison-linear 
equations for 18 unknowns. 

It should be noted that successive orientations should be as different 

1 2 

as possible. For example if R and R are very close to each other, very 
little new information would be provided by a second set of measurements. 

In the presence of errors it would then be impossible to solve the equations 
for (/() and (»}). 


If we expand the ij^^ element of equation (35) we obtain 




-1 


U=1 




xj 


(36) 


where a - 1, 2, ...N and i and j = 1, 2, 3. 
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In a subsequent section we show how this set of equations may be solved by 
an iterative procedure using a Taylor series expansion which linearizes 
equatiwi (36). 


A procedure for determining the dii~ection cosines of a single 
calibration coil and one sensor 


The solution of the set of equations defined by equation (36) is a 
ct^lex procedure. Furthemore it requires that the orthogcwial sensor array 
be rotated relative to the calibration coils. In some cases this may not be 
possible to do. Hoi««ver, if we can determine the direction cosines of the 
calibration coils in separate experinents we can use the procedure described 
above to obtain the magnetic axes of the sensor. In this section we 

will show how a single sensor may be used io detemins the direction ccsines 
of one calibration coil. In subsequent sections we generalize this method to 


the case of three sensors and three coils. 

Let n. be the unit vector defining the magnetic axis of the sensor in 

yv 

sensor coordinates, n-the unit vector defining the calibration coil axis ’*n 
geographic coordinates and B- the magnitude of the calibration field. The 

«l 

output of the sensor is given by 


M. . = tn. • B. 

ij 1 

where B. must be defined in sensor coordinates. But B. = B.n.' with n.' 
-3 J J J 

being the representation of n. in sensor coordinates. Thus 

tJ 


Hij/Bj . n,. . n.' 


A ^ 

considering m- and n.' as coluirn vectors we can rewrite the dot procuct 

* J 

M. ./B. = ^ m . n . » = £ (m^). . n 


ij j 




zt. ot 





43 


But V(SEN) = (R) y(GEO) 
so that 


or 


n*. = ^R.n-- 

f 


Hence 


or 


M. ./B. = i i{m). R^.n. . 


R^o "5- = M../B. 

tfi p kt Pj ij' j 


(37) 


This equation is identical to eq. (36) when the calibration matrix (B^) 
is diagonal, i.e. only one calibration coil is excited in any one measuranent. 

Since i and j are fixed there are six unknowns in eq. (37). To solve 
for them ;« inust perform at least six different experiments with the sensor 

A ^ 

m- taking different orientations relative to the calibration coil n.. Using 
(0 to designate each of these experi.T.ents we have 


^i "(3 = «i3<^^'^®3^^^ 


)T = 1. 2. ... N 


(38) 


To linearize this equation we assume that we know the orientations of 

and n. to about one degree. This is reasonable since both the coil axes 
3 

and magnetic axes can be manufactured with this accuracy. Thus we take 
initially. 


m . 
cn 


m 




. o 

f) ^ n *t • 

P3 PJ 


(39) 


We then make a Taylor series expansion of the left hand side of (33) about 
this initial model, obtaining 

flX;, (40) 


IJ J A-I 


'>0 

X 
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In this expansion x is a six c(x;^onent column vector constructed frtxs 
the two unknown vectors, i.e. 


Thus 


Hence, Aj(?>) » n?j 
The derivatives are given by 


- & 
.0 _o 


3 : the initial cK)del 
- corrections to the initial model 


55^ 


A 

Xo 


= if IcsO) 

•< f *”Ai *•> 


*0 

X 


but 


"%• . ^ 
J5^ A 


so that 


- i "°5j 


AO f 
X 


Similarly we find 
3A 




X 


Substituting into eq. (40) gives 


Mj ,(?)/8 (5) - f £ R 5(.) n®3 ♦ £ (£ R n®. ,ym 

•J J at 5 • £Xi v’l ;^i.| ^ 41 




A\ 

We may simplify eq. (41) by the use of the matrix notation. 


( 41 ) 


»ij(« “ 


Let 
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Then 

bjj{J) - (S®j)^.(R(I)ii/) ♦ (R(8)^/)-iiS, ♦ (R^(?)nij“)-/lnj («) 

This equation nay be interpreted physically as follows. The left side 
is the norralized Rieasurement of the i^^ sensor in response to the cali- 

bration coil. The first term on the right hand side is the predicted, 
normalized response when we use the initial model for the orientations of 
the sensor and coil. The second tens on the right hand side is the projection 
of the initial coil axis on the correction to the sensor axis in sensor coordi- 
nates. The third term is the projection of the initial sensor axis on the 
correction to the coil axis in geographic coordinates. 

If we iiX3ve the first term of the rightside to the left we obtain 

{R(})S®MS. ♦ (RT(i!)m'>).(!Sj . b.j(0 - t(£=r.RU)A“] for S'" 1 . 2. ... ti (43) 


In this fom, we clearly have a set of N linear equations in the six unknowns. 


4m, 4n. The fonn of this equation is 

3 

£ X. = Yv If = 1, 2. ... N 


(44) 


where are the elonents of the correction vectors 4rn, fin; is the difference 
between the measured and predicted values of the normalized measurements; and 
Ajp, is an (N x 6) matrix with rows constructed from the initial coil axis 
in sensor coordinates and the initial sensor axis in geographic coordinates. 

We solve the set of N eq's (43) in an iterative manner. An initial guess 
is made for the direction cosines of pI| and nj. By construction these are 
designed to be as nearly aligned v/ith the corresponding axes of sensor and 
geographic coordinates. Thus the initial guess is 

(Sj“) =S. and (A“) 

Evaluating the constant terms in eq. (43) we obtain a set of equations of 
the form (44). These equations are solved for the correction vectors fim ^ 
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1 

fln • by the method discussed below. The corrections are added to the initial 
guess to obtain a new starting model. Thus 


*i * ® I 
AO 

n. ® n I + An . 


Using these new vectors as a starting model, we repeat the procedure obtaining 
a second set of solutions m- , n/. Provided the initial model is close to 

* «l 

the correct model, the problem is nearly linear and the sequence of corrections 
converge rapidly to zero. 
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Solution of the set of H equations and M unkno’.<ns 

In the preceding section we showed the problem of simultaneously 
determining the direction cosines of sensor and coil could be reduced to a 
problem of solving N linear equations for six unknowns. Thus we solve 

Ajp H 3f= 1. 2, ... N (44) 

for the m (m = 6) unknowns X^. In matrix form this may be written 
(A)x = y 

A general procedure for doing this has been described by Lanazos, 1961 (p. 100- 
162). We briefly summarize this procedure below. 

Let (A) be the (N x M) matrix of coefficients in equation (44). Construct 
the (N X H) syiiaretric matrix AA^ and tne (M x M) syiniietric matrix A^A. Solve 
the two eigenvalue problems 

AA u ='T u N eigenvalues 

T kl 

A A V = ^ V M eigenvalues 


Construct two matrices U and V from the eigenvectors u and v. The matrix 
U is (N X N) and V is (H x M). At most there will be only p, non-zero 
eigenvalues <1- where p< min (N,M). The non-zero eigenvalues are the sane 
for both AA^ and A^A 

Now construct an (N x P) matrix Up and an (N x (N-P)) matrix by 
partitioning U into two matrices having columns corresponding to non-zero 
and zero eigenvalues respectively. Similarly construct an (M x P) matrix 
Vp and an (M x (M-P)) matrix from V. The matrices Up, Vp diagonalize 
the matrix A by the transformation 



(46) 


where A_ 
P 


is a p X p matrix with diagonal 


elements corresponding to the square 
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roots of the non-zero eigenvalues of eq. (45). From (46) we have 

A=Up.ApVp^ (HxM) (47) 

The natural inverse of (A) is given by 

B = VpAp"’ UpT (M X N) {48) 


AB = UU^ 

BA = VV^ 

The solution to our set of equations (44) is then as follows 


(49) 


Ax = y 

(UpAp Vp^) X - y 

* ‘ <Vp'’ 

X = (B)y (50) 

The solution, eq. (50), always exits, but it will not necessarily be 
a good one. First, sane eigenvalues may be non-zero but snail. In this 
case, corresponding eleinents in the inverse matrix -A^p”^ will be very large. 

As a consequence, very small errors in the constant vector y will be magnified 
in the solution vector x. 

A second problem is that errors made in the measurerrent of A and the 
vector y may make the equations incompatible. Lanczos shows the set of 
equations are compatible only if 

y ' 0 (51) 

Physically this implies that the constant vector y must be orthogonal to 
all eigenvectors corresponding to zero eigenvalues. 

Finally, the set of equations may be deficient. In this case they contain 
insufficient information to determine certain linear combinations of the 
unknowns. These combinations are given by 


( 52 ) 
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where is the matrix of zero eigenvectors from the matrix V and is an 
arbitrary column vector. 
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Co mputer si m ulation of the simultane ou s de termi nation of the 
orientations of a single sensor and singl e coil - tio 
constraints 

To test the procedure described in the preceding two sections, wa have 
used a computer program to simulate an actual calibration experiment. Inputs 
to this simulation program are the direction cosines of the coil and sensor 
in their respective Cartesian coordinate systems and six transformation matrices 
corresponding to different orientations of the sensor relative to the coil. 
Outputs from the simulation program are the normalized magnetometer measure- 
ments which would have been made in the absence of measurement errors. 

Six different sensor orientations were used as summarized in Figure 8. 

For this experiment, the sensor was chosen to be along the x axis of sensor 
coordinates and the coil along the z axis of geographic coordinates. The 
nine different transformation matrices and the corresponding normalized sensor 
mieasurements are presented in Table 1. A listing of the simulation program 
is included as Appendix A2. 

The second step in the simulation required a computer program which 
implements the procedure described in the preceding two sections. This 
program was written in the IBM-TSO (Time Sharing Option) version of SPEAKEZ. 

This language v/as written and is maintained by the Argonne National Laboratory 
(Cohen and Pieper. 1976). It is especially designed to facilitate the manipu- 
lation of vectors and matrices. A listing of the latest version of this 
program (MAGCALl)is included as Appendix A3. 

The results obtained by this program were quite satisfactory, converging 
to nearly correct values after at most three iterations, A representative 
result using five measurements 

= (.99935, .030011, .039929) 
as compared to the kno'.vn input, 

M - (.99375, .030000, .040000) 

X 
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While the foregoing procedure is successful it is difficult to carry 
out experimentally. The problem is that it is exceedingly difficult to 
perform the measurements necessary to calculate the transformation matrix 
(R). Thus tne most desirable procedure is one which minimizes the number 
of different orientations of sensor and coil. 

One method of reducing the number of necessary orientations is to take 
advantage of the unit vector constraints. Since both m. and n. are unit 
vectors there are only four unknowns rather than six. Hence only four orienta- 
tions should be required. 


Utilization of the unit vector constraint to reduce the 
number of experimental measurements 


Thus far we have not used the constraint that both m. and n. should be 

• w 

unit vectors. This constraint can be included in two ways. The most straight 
forward v/ay is to linearize the constraint equations and include them in the 
set of linearized equations solved by the computer procedure. A second way 
is to use the constrain^" directly to eliminate two unknowns from eq. (38) 
and then linearize these modified equaf-ions. We find empirically that only 
the second method works. Basically, the reason appears to be a result of 
losing too much information in separately linearizing the fundamental 
equation (33) and the constraints. 


To include the unit vector constraint directly we return to the original 
equation describing the result of any particular measurement, eq. (38). 

ir' 1. ... N 

Since iri and n are unit vectors we have 


f(m .)^ = 1 and ^ = 1 

^ --V ""J 


( 33 ) 



5Z 


we may solve these for the diagonal elements 


"■ii - X 


'Vi>' 


n,, = Jl - 


(53) 


(n .) 

while these can be substituted directly in eq. (38) it is not necessary to 
do so. Instead, we treat eq. (38) as before as a function of six variables, 
i.e. 

f (m n.,n.,n.) 

' XT yi’ zi’ xj yj* zj' 

but utilize the rules for implicit differentiation to evaluate the derivatives 
of the diagonal elements. Thus, 


f Hi.) * 


where 


df = ? dm . + f dn . 


Isolating the diagonal element in each sum gives 


df 




dm . + 






— dm . . + £ 

L ■ n ^ . 
n oc- « 


?f.d„ .^JL 

2n.. oc.i tfn.. 


. <? Mj • 


dn 


JJ 


jj 


(54) 


(55) 


(56) 


From equations (54 & 55) the total derivatives of the diagonal elements are 


^ 3m.. 

dm. . = ^ -5—^“ dm . 

n 3m^,- 

/ 3n . . 

dn . . - ^ dn 

W 0(1;) ^ 

Substituting in eq. (56) gives 


df = ^.( 


-v - 1 ‘^m ,• 
‘ oil 


2m- ■ 2m cO 


£. (v 


Pf 


ii "’oti 

.Pf 


' '"jj 


P n . . 

^— ^) dn . 


</n 


(57) 




Since takes on only tv;o values in each sensor we have clearly reduced the 
number of unknowns from six to four. The remaining two unknov-/ns are defined 
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by the equation (53). Since 


f- 

we may evaluate the derivatives in (57). Thus 

^7 “ j, "oU 


!!!ii 

?ro . 
>1 


% 

"ii 


^ "jj 


c^i i 




Hence the quantities within parenthesis in equation (57) are 
?m.. <?m^. 3cC J ' '"ii 


1 ^ + ^, ^3 = C(R^m,-L - (R^m-). (^)3 
vn ■ <?n.. on^.-* i be ' I'l 'n..''-' 

<?cJ JJ 


n . 

SA; 

Jj 


( 53 ) 


Tpese derivatives must be evaluated for the initial model, i.e. m. and n- 

^ 3 

are replaced by the guesses m.° and r .° in the above expressions. 

* w 

Then, proceeding as we did in the derivation of eq, (<13) we find 


s<f t 


"’ii 


O 

n . 


t - (R^m °) (^^)] fin 


JJ 


' Mjj/Bj - rS°] 


(59) 


Note for convenience we have not indicated the dependence of R, M.. and B- 

1 J J 

on (j), the relative orientation of sensor i and coil j. 

This exp>“es 5 ion is very similar to eq. (43) except it no longer has 
an obvious physical interpretation because we have explicitly included 
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the unit vector constraint. The major difference is the reduction of the 
number of unknowns from six to four. Thus the equation now has the form 

(A) X = y = 1, 2, ... N fcO ) 

These equations may be solved by exactly the same method as used before. 
The only changes required are in the dimensionality of the (A) matrix and 
the calculation of its elements. 


Ccxsputer simulation including constraints 


The procedure described in the previous section was implemented by 
making slight modifications in the program FIAGCAL. A listing of the new 
program MAGCiU.1 is included as Appendix A3. This program was tested 
using the same data as were used to test KAGCAL (c.f. Table 1). On a basis 
of this test we conclude the second procedure is far superior to the first. 

He find that the inclusion of the unit vector constraint improves the 
rate of convergence, inproves the accuracy and reduces the number of necessary 
measurejnents. A san^le result obtained after two i.erations was 
= (.99875, .030000, .039996) 
n^ - (.010001, .020006, .99975) 

These should be coc^ared to the known values 
= (.99875, .030000, .040000) 
n^ = (.010000, .020000, .999750) 

Similarly good results were obtained using only four experiments, provided 
experiment number 3 of Table 1 was not included. This experiment corresponds 
to alignment of the sensor with the calibration coil. 

This UP' xpected result suggests that some sensor orientations are better 
than others :or determining the direction cosines of sensor and coil. 

Ke interpret this result in the following fashion. The function 


f - 


nip 


m 


«■« 
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A ^ 

is a naxirnum for a matrix R corresponding to alignment of m. and n.. Our 

R£thod of linearizing the equation to solve for m-, n. utilizes the deriva- 

* J 

A A 

tives of f as a function of the elements of m and n to calculate corrections 
in the initial guess. Since cur initial guess corresponds to exact alignment 
the inl».:al de'-ivaLives should be exactly zero. Thus there is no way to 
calculate a correction to the initial guess. 

Examination of the eigenvalues obtained in the solution of the set of 
four equations for this cawuter siniulation shows that one of the four eigen- 
values was much smaller than the remaining three. This implies that there 
are fewer equations than unkncvins and therefore the equations are not really 
soluble. Clearly the data obtained in experiment 3 by near alignment of 
sensor and coil does not contribute much information. 

Errors in magnetic axes orientation using simultaneous determination 
procedure 

In our procedure for simulatneous determination of sensor and coil axis 
we utilized the Lanezos inverse. Thus 
Ax = y 

has a solution 
x = By 
where B = 

As discussed by Jackson, 1972, the errors in the solution vector x are given by 
Var (x^) = i Var (y.) (61) 

Since depends on the reciprocal eigenvalues of Ap * any small eigenvalues 
will cause very large errors in all components of the solution vector. 

Equation (61) may be written in matrix form as 
var (x) = (B^) Var (J) 


( 62 ) 
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2 

In this relation B is a matrix whose elements are the squares of the elements 
of B- Var (x) and Var (y) are vectors constructed from the respective 
consonant variances. 

To obtain some idea of the errors to be expected in our simultaneous 
determination procedure we have included eq. (62) as part of program MAGCAL 
and MAGCALl For the simulation described in the preceding section (four 
experiments using constraints) we find the following results. 

First we assume the variance of all input measurements to be the same, 
and roughly of order 1 garana in 50,000 /, i.e. VAR(y^ (1/50,000) . 

Then from eq. (61) 

Var (x^) = VAR(y) • £ (B^-)^ 

2 

A typical row in the matrix (S ) was for our first simulation 
(.5. 3.0, 1.0, 5.5) 

which sums to 10. Thus the expected variance in the solutions is roughly 

ten times that of the input data. Hence 

VAR (x^) 10 (1/50,000)^ 

-5 

Or RKS error in X|^ -^6 x 10 . 

To convert this to angular measure we note 
cos ^ = n^^ = .02 i 6 X 10~^ 
and 

^0 - Arc cos(.02 +6 x 10~^) - Arc cos (.02) 
or 

^ 3.4 X lO”"' = 12 arc seconds. 

This error is quite large compared to the design goal of one arc second. 

One possible means of reducing this error is to reduce the magnitude of 
the elements of the inverse matrix. Since this depends on the original 
matrix (A) which in turn depends on the relative orientation of sensor and 
coil, it is possible that a different set of four experiments would provide 
a more accurate determination. 
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To check this possibility we have performed a second simulation using 
the four orientations shown in figure 9 and Table 2. To decide on which 
orientations to use we guessed that large elements of the inverse B depend 
on small elem'-nts of A. Apart from zero elements, the smallest elements of 
A arise from the 45® relative orientations. We thus decided to use only SO® 
rotations to obtain the four orientations shoivn in figure 9. 

The results from this second simulation were quite surprising. We 

found 

M AX (RMS (x,^)) - 2.5 x lO"^ 

This is more than a factor of two improvement over the first set of four 
orientations. The corresponding angular accuracy is of order 6 arc seconds. 

While an error of six arc seconds does not meet the design goal it 
is at least as good as can be expected if the coil cosines are determined 
with a conventional magnetic theodolite. This error might be fu. ther reduced 
by using the results of an initial determination to design an optimum set of 
four measurements. It could be reduced still further by repeating the set 
of four measurements a number of times. 
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Calibration of Three Sensors Using Three Calibration Coils 

In the preceding section we shaded how a single sensor could 
be calibrated using only one coil. This procedure requires a sufficient 
number of relative orientations of the sensor and coil to define all 
of the unknowns. In this section we show how the use of three coils makes it 
possible to calibrate three sensors with the same number of orientations. 

Let us assume as before, that for each sensor orientation we carry out 
three measureitients. Each measurement is the vector output of the sensor 
assembly for a given calibration field. The expected fields for this set of 
three measurenents may be written 

(B„) = R n] (BJ (3C) 

m V, 

where B and are 3x3 matrices with column corresponding to the three 
m c 

measured and calibration fields respectively This may be written 

Rn = b (63) 

where b = (Bjj^)(B^)~* is the "nonnalized" observation matrix. 

To linearize this set of equations for p and n we assume an approximate 
model is knew, i.e. y = y^ and n = We then v/rite y and n as small pertur- 
bations about this known model, thus 

u = Pq + 5y 
n = Hq + 5fi 

Substituting in equation (63) we find 

(u^^R)<in + ^^^0^ = [b - p^^ Rn^] (64) 

where we drop the second order tern 6y^6n which we assume to be very small. 
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The righthand side of equation (64) is the difference between the actual 
measurerrents (b) and the measur«r.ents that would be expected as a result of 
the initial nocel- The unknowns 5n. on the left hand side are the 
corrections to the initial model required to obtain better agreement between 
the observations and the model. 

This set of nine equations in eighteen unknowns has the form 

(A) X = y (65) 

Here x is a colunm vector with 18 rows constructed from the unknown elements 
of 6y and 6n. The vector y is a column vector with 9 rows constructed from 
the residuals between the measured and predicted magnetic field measurements. 

The matrix (A) is a 9 x 18 rectangular matrix with elements constructed 
fran the two matrices (uq^R) aod (ROq). Since this matrix equation consti- 
tutes only 9 relations betv;een the 18 unknowns we must obtain additional, 
linearly independent relationships. We do this by using a different orienta- 
tion of sensors and coil, l-rovided these orientations are properly chosen 
the coefficient matrix will be non-singular and the set of equations will 
be soluble. 

To demonstrate the feasibility of the above app»*Odch we have carried 
out a conputer simulation as was done for the single sensor, single coil 
calibration procedure. Direction cosines of the magnitude expected for sensor and 
coil were chosen arbitrarily (Table 2). Then equation (63) v/as used to 
calculate the expected magnetic field measurements. Seven different 
orientations of sensors and coils were used as this v/as the minimum necessary 
to define all sensor axes wiien only one calibration coil is used. The seven 
transformation matrices and the seven calculated measurement matrices v/ere 
then used as data in the perturbation procedure described above. 
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To set up the set of nine algebraic equations corresponding to the matrix 
equation (64) v« take the ij^^ elercnt. 



Clearly, the equation couples only the coil axis with the i^^ sensor, 
i.e. each equation links only six of the 18 unknowns. Thus if we adopt the 
convention of first fixing j and then allowing i to run through its range 
we generate the set of 9 equations shown schematically in Table 5. In 
Table 5 each column contains the coefficients of the three unknowns associated 
with a given unit vector. Each row corresponds to successive values of j and 
i. The entries within this table, i.e. ROW 1, COL 3 refer to rows or columns 
of the model matrices, (y„^R) and (Rn ) shown at the bottom left and right 
sides of the Table. All blank entries correspond to coefficients of zero. 

This representation of the set of nine eqjations is particularly convenient 
for the computer language, SPEAKEZ. Because this language allows the manipu- 
lation of vectors and matrices as entities it is possibly to efficiently 
define the 9 x 18 coefficient matrix (A) simply by placing the appropriate 
rows and columns of the model matrices at appropriate locations of (A). 

A program v/hich generates the (7*9 x 18) "A" matrix corresponding to 
the seven orientations of sensor and coil sho'.^n in Table 3 was written. 

The solubility of this set of 63 equations was then tested in the following 
manner. Assume an initial model corresponding to exact alignment of both 
sensors and coils in their respective coordinate systems, i.e. ^^and are 
identity matrices. Furthermore assume the transformation matrices correspond 
to exact alignme.its of the sensor and geographic coordinate systems (not 
normally the case fo»- actual measurements). With these assumptions the 
columns of (ROq) sre the unit vectors of the coils expressed in sensor 
coordinates. However, because we assume the coils are exactly aligned with 
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geographic coordinates the columns of (Rn^) are actually the unit vectors of 
geographic coordinates as seen in sensor coordinates. Similarly, rows of 
(y^^R) correspond to columns of (R^p^), but by the same arguments as above 
these are the unit vectors of sensor coordinates as seen in geographic 
coort. lates. Using these facts we can inwiediately write down the elements 
of the A matrix by inspection. To decide whether this set of 

equations can be solved we take the determinant. If this is non-zero the 
matrix is non-singular and the equations soluble. 

For the A matrix based on the seven orientations of Table 3 SPEAKEZ 
methods showed the determinant was zero. Thus the equations were insoluble. 
Our interpretation of this result is as follows. The seven orientations of 
table 3 vere deliberately chosen to produce a specific set of four orienta- 
tions of each sensor relative to one coil. These sets of four orientations 
were chosen so that the four unknowns associated with one coil and one sensor 
could be determined. Only four unknowns were present because the unit vector 
constraints ivere used to eliminate two of the six unknowns. 

This argument leads us to the conclusion that at least several 
additional orientations would be required to solve for the 18 unkncn<ns by 
the above method. Because of the experimental difficulties associated with 
the measurement of the R matrix we rule this out as a viable procedure. 
Consequently, we extend the foregoing procedure to make use of the six unit 
vector constraints. 

The constraint equations may be written 



I n / = 1 

a=l 

3 p 
I tn = 1 
a=1 


( 67 ) 


Solving for the diagonal elements v/e have for columns of (n) 
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n.. = 
JJ 


1 - J n 
aj*j 


ay 


In our perturbation approximation 


n . = n . + fin . 
aj aj oj 


So that the diagonal elements become 


n. . + 


in,, = [1 - r (n + in„,)‘] 


2i1/2 


JJ ■ J3 " ' aj ’ aj' 

Expanding the square and dropping second order terms gives 


n,,”+in„ [( 1 - I. (n„,V )-2 I n in .] 


i1/2 


JJ 


JJ 


O^j 


aj 


aj'j 


aj 


aJ 


( 68 ) 


But 


Ox2,I/2 


-jj ■ ” - J > ' 


so that Q 

"i/ ^ = "j/ 1’ - ^ 

Because the second term within the brackets is much smaller than one we 
have 


y n . fin 

fs . I 

af3 


aj aj 


n.. + fin.. * n.. [1 - 

JJ JJ JJ *■ (n..°)^ 


-] 


JJ 


or finally 


n . 
oJ_ 


fin.. ~ - I {-^) «n . 


(69) 


In an analogous fashion 

0 


fim 


m . 

• - = - I (-^) 

n j, _ o I 


a/'j m. 


a1 


11 


(70) 
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To include these constraints in a modified A matrix we note the following. 
The two constraints may be written in the form 



Now multiply the first equation by the colunai of the A matrix corresponding 
to the unknown 5n... Similarly multiply the second equation by the column 

wO 

of A corresponding to 6m^-. Then, add the two equations obtaining 
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Finally, add this equation the ij^^ equation previously obtained without 
using constraints, (eq. (66)). 

When this procedure is carried out explicitly it can be seen that it is 
exactly equivalent to substituting the constraint equations for the diagonal 
elements and then rearranging terms. Thus to modify the A matrix shown in 
Table 5 we must construct a (9 x 18) N matrix with the elements shown in 
Table 6. Then this matrix (N) must be subtracted from the matrix (A). 
Subtraction will introduce zeros in the columns of the modified A matrix 
corresponding to the diagonal elements of the unknown matrices (5^), (5n)- 
If we now eliminate all comumns of the modified A matrix corre.»ponding to 
these diagonal elements we will obtain a new set of equations of the form 

(AA) x' = y 

where now x' is a column vector with 12 unknowns taken from the off diagonal 
elements of (<5n), (6u). y is a 9 row column vector of residuals between 
observations and predictions as previously defined. The matrix (AA) is a 
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(9 X 12) rectangular matrix obtained as described above. 

A SPEAKEZ program "EQUATS" was written to implement this procedure. 

A listing of this program is included as Appendix A5. in the initial text of 
this program all seven orientations of table 3 were used to construct 
a (63 X 12) matrix (AA). The determinant of this matrix was found to be 
non-zero for the particular simulation used. Consequently two additional 
programs were written to solve for the direction cosines in an iterative 
fashion. 

Program "SOLEQU" (listing also included in Appendix A5), utilizes a 
SPEAKEZ linkule (similar to subroutine) to solve the set of equations defined 
by the matrix (AA) and the vector y. This linkule utilizes the singular 
value decomposition procedure of Lanzcos described earlier. 

Because this linkule is w/ritten to minimize time and storage requirements 
it does not provide as much information as our own implementation of this 
procedure embodied in program MAGCAL and MAGCAL 1 (See Appendix A3). 

Ho'wever, because our implementation utilizes too much storage we were unable 
to use it to solve the set of 63 equations. 

A second program "MAINPO" was written to implement the iterative perturba- 
tion procedure our method of simultaneous solution is based on. This program 

initiatizes the calculation first reading in the measurements and transformatio 

0 o 

matrices, and then reading a first guess of the unknowns, (n ) and (p ). 
Normally we assume these are identity matrices since these are the design 
goals for both sensors and coils. The program then calls program EQ'JATS. 

As its first step program EQUATS increments (y°) and (n°) using previously 
calculated corrections (6y) and (^n). (Initially these are zero). EQUATS 
then uses the constraint conditions to calculate the diagonal elements of the 
modified matrices (w° + 5p) and (n° + <5n). It then proceeds to generate 
the matrix (AA) and the vector y. Control is passed to program MAINPO which 
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calls SOLEQU to solve the set of equations for the vector x. In its final 
step SOLEQU creates the correction matrices (6ii) and (6n) from elements of the 
solution vector x. Control •“ returned to program MA INFO which loops back 
to call program EQUATS. This procedure is repeated the number of times 
specific .during the initiatization phase of program MAINPO. 

Results obtained with program HAIMPO are comparable to those obtained 
with the programs which were written to solve the case of one sensor and one 
coil. For example, the second corrected model was equal to the model used in 
the simulation program to better than the fifth decimal place. 

To determine the minimum number of orientations required to carry out 
a procedure we progressively reduced the number of different orientations 
input to program MAINPO. When this nu.iber was less than four, the determinant 
became zero and we were unable to obtain a correct solution to the proolem. 
Although we have only examined the case corresponding to the first four 

■-‘entations of table 3, we believe that any four orientations satisfying the 
foil: wing criteria v/ould be sufficient to define the 18 direction cosines of (u) 
and (n) using our procedure. These criteria are the same as was found for 
calibrating a single sensor with a single coil: (1) there must be 4 different 
orientations of the sensor and the coil, (2) direct alignment of the sensor and 
coil provides no useful information, (3) rotations of 90° from the initial 
position generate satisfactory orientation provided one rotation brings the sensor 
axis initially transverst to both coil and sensor parallel to calibration coil. 

In surmiary, our procedure for calibrating the direction cosines of a 

three axis sensor assembly using a three axis Helmholtz coil system is the 
following. Attach an optical cube to the sensor assembly. Place the sensor 
assembly in a fixture that allows 90° rotations about two orthogonal axes 
and which does not allow the center of the cube to translate during rotation. 

Begin with an alignment which places the x', y', z’ axes of the sensor, i.e. 
of the optical cube, in near coincidence with the x, y, z axes of the coil, 
i.e. of geographic coordinates (c.f. Figure 9 )• Level the fixture using 
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two theodolites in the horizontal plane on the positive x and y axes of 
geographyc coordinates. The leveling procedure should be sufficiently 
accurate to guarantee that theodolite cross hairs remain reflected within 
their fields of view after successive 90° rotations about the two axes of 
the fixture. Use the theodolite setting circles to measure the actual 
orientations of the x' and y' axes of the optical cube to one arc second or 
better. Apply an accurately known field in the x calibration coil and measure 
the output of the three sensors. Repeat for the y calibration coil and the 
2 coil using exactly the same field magnitude. From the measurements 
construct two matrices (8^^^) and (B^) using the successive measurements and 
calibration fields as columns. Thus. 




Next, calculate the normalized measurement matrix 




which by the choice of calibration fields becomes 



Now create a nine element column vector of normalized measurements by placing 
the three columns of the matrix (b) successively beneath each other. Thus 
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Finally, using the orientation angles determined by the theodolites construct 
the transformation matrix (R) from geographic to sensor coordinates. 

The (3 x-3) rotation niatric (R) and the nine element column vector 
M are the data to be input to the simultaneous equation solving program MAIIIPO. 

Having completed the measurements for the first orientation we rotate 
the sensor assembly + 90° about its z axis (azimuth axis of fixture). The 
actual orientations of two sensor axes of the sensor assembly are determined 
by the x and y theodolites as before. We then repeat the previous sequence 
of three calibration fields and corresponding measurements. A matrix 
and a vector M calculated as above form the second set of input data * / u- , 
MAINPO. 

For the third orientation of figure 9 we rotate an additional + 90° 
about the azimuth axis of the fixture and obtain rneasurements as before. 

Finally, the fourth orientation is obtained by a 90° rotation about the 
elevation axis of the fixtur'e. This brings the top of the optical cube dov,r. 
to a position where it may be viewed by the y theodolite. 

After all data have been properly entered into program MAINPO, the program 
requests an initial guess for the tritrices (p) and (n) defining the direction 
cosines of the sensors in sensor coordinates and of the coils in geographic 
coordinates. If no prior information is available these are assumed to be 
identity matrices. The program then proceeds through several iterations 
obtaining a final, least square determination of the best fit of the direction 
cosines of sensors and coils to the input data. 

Recommendations 

The high absolute accuracy desired for the MAG5AT magnetometer requires 
accurate calibration equipment and some new procedures. We recortf'^nd acquisition 
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To make an accurate determination of sensor direction cosines, the 
calibration coil direction cosines must also be determined accurately. This 
can be done in tivo ways. 

1} With a magnetic theodolite which would need to be purchased 

2) Sirnultaneously as the sensor cosines are determined 

In our report we choose the simultaneous determination procedure. We 
show that this technique should be able to provide accuracy of order 5-10 arc 
seconds. This accuracy is comparable to that of easily acquired magnetic 
theodolites. To significantly improve this accuracy would require repeated 
measurement of the direction cosines. This procedure requires four reorien- 
tations of the sensor assend>ly. For each orientation both theodolites imist 
be read and a series of three calibration fields applied and the magnetometer 
rvtput recorded. 

We recommend that this procedure be automated and placed under the 
control of a minicomputer. This would require that the rotation fixture have 
motor control of each axis. Additionally, the theodolites should have digital 
readouts. 

If we accept 10 arc seconds as the obtainable accuracy of the direction 
cosines of the sensors, then the absolute accuracy of the magnetoneter cannot 
be better than + 2 gamma. However, individual sensor errors will also be 
about + 1 gamma, so we expect the final accuracy to be about + 3 gamna. 
Significantly improving this accuracy appears to involve very labor intensive 
calibration procedures. 
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of the following equipment: 

1) Two optical theodolites of one arc second accuracy 

2) Six digit digital voltmeter 

3) Two precision levels of one arc second accruacy 

4) Snx)oth. non-magnetic tost table of granite or glass 

5) Brass fixture for rotating sensor assembly about horizontal and 
vertical axes with precision of order one minute 

6) Minicomputer controller and data logger with video display keyboard 
entry 

Several ax>difications of the existing facility also appear to be 
required. These include: 

1) Separate pillars sunk to bedrock on the North and East calibration 
coil axes. These will hold the two theodolites. 

2) Two optically flat ports at right angles to each other in the thermal 
shroud to monitor changes in sensor orientation as a function of 
temperature. 

One raajor problem likely to be encountered ii? the calibration procedures 
is digitizing the magnetoneter output with sufficient precisioj. If a 16 bit 
Analog to Digital Converter is used for the +64K to -64K gamma dynamic range 
the least significant bit will correspond to 2 gamma. This uncertainty 
of + 1 gamma limits the accuracy of the calibrations to an unacceptaboe level. 
Two solutions for this problem are possible: 

1) Use a data logger to make a large number of samples and average 

to higher precision, (This works only if the facility noise is large.) 

2) Use AC signals rather than DC signals as the calibration fields. 

If facility noise is low, small signals must be present in calibration 
axes orthovjonal to main calibration fields. Use least square sine wave 
fits to determine AC signal amplitudes precisely. 
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Figure Captions 

1. Schanatic illustration showing the ar rangeunent of calibration coils 

test table and theodolites. Theodolites and levels define the orthogonal, 

. . V V C The direction cosines of the calibration 

geographic coordinates X, Y, Z. 

coils in geographic coordinates are given by the non -orthogonal unit 

A ^ A 

vectors n^. n^. 

2. Diagram showing a possible arrangement of three ring core sensors 
supported by a plastic cube and topped by a silvered, optical octagon. 
Normals to the faces of the octagon establish an orthogonal sensor 
coordinate system. The direction cosines of the nagnetic axes of the 
ring cores in sensor coordinates are given by the non -orthogonal unit 

A A A 

vectors 13^^, iiiy, in^. 


3- A photograph of an earth inductor (or magnetic theodolite) showing 
the type of fixture required for calibrating the angular orientation of 
a sensor magnetic axis in sensor coordinates. As discussed in text the 
rotating search coil is replaced by a mounting plate to hold sensor. 


4. A schematic diagram showing the modi f ications of the earth inductor 
required to obtain a magnetometer test fixture. 

5- A schematic illustration showing how angles 0 'y 90° and ^ 0° 

are measured by theodolites to determine absolute orientation of normals 

to the faces of an optical octagon. 


6. A fixture for aligning individual sensors with the calibration field. 

7. Typical noise spectrum for a ring core fluxgate sensor. Horizontol 
brackets at top show frequency band measured by an experiment of given 

duration and sample rate. Dashed horizontol lines show qu.Tnt i zat ion 
noise for cB = ly and 53 % l/32y. 



8 . Six relative orientations of sensor and coil used in first 
computer simulation. 

9. Four relative orientations of sensor and coil necessary for 
Full calibration of three sensors and three coils. 
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First Computer Simulation 
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List of Table s 

I. Transformation natricies used in simulated calibration of a single 
sensor with a single coil. 

T T 

7., Direction cosines of sensors (u ) and coils (n ) assumed in 
simulation of full calibration. 

3 . Transformation r -t icies for seven different relative orientations 
of sensors and coils in full calibration. 

h. Normalized creasurements for seven different orientations of sensors 
and coils. Reading across and down the cieasurements correspond to 
outputs of sensors I, 2, 3 for a field first in coil 1, then coil 2 
and finally coil 3- They then repeat for the next relative 
orientat ions. 

5 . Coefficients in the set of 9 equations produced by a fixed orienta- 
tion of sensors and coils. Rov< I and Col. J refer to the two 
Riatricies shown at the bottom. 

6. Modifications to be subtracted from the coefficients of Table 5 
caused by the use of constraints. The notation A 318 refers to the 
9th rov; and l8th column of Table 5, and so on. 
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Al. MAGSAT oagnetOMter specifications. 

A2. SPE/UCEZ program "NAGSiif' which generates siowtated sensor output 
when excited by a calibration coil in a given orientation. Inputs 
are the direction cosines of the sensor and coil in their respective 
coordinate systems. 

A3* SPEAKEZ progr^ "HAGCAL 1*' which solves for the directicwt <x>sine$ 
of sensor and coil using initial guess and simulated s^sur^ients. 

A4. SPEAKEZ program "MEUSIH" which generates sirajlated measurements 
for three sensor and three coil simulation. 

A5* SPEAKEZ programs •‘KAIMPO,” "EQUATS,** and “SOLEQU." Together these 
programs use measurements of the transformation matrix and magnet- 
ometer output with an initial guess to solve for direction cosines 
of three sensors and three coils. 
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Appendix A1 

MAGSAT magnetaneter specifications . 

Specifications for the Measurement of Vector Magnetic Fields 

1. General - The magnetometer system described in this specification shall 
^asure three orthogonal exponents of the ambient vector magnetic field aboard 
a l<w altitude earth-orbiting spacecraft. It is specifically intended for 
flight on the NASA Applications Explorer Mission spacecraft. The 

spacecraft Mill be launched by a Scout vehicle into a near polar, low altitude 
orbit for the purpose of making a global survey of the earth vector autgnetic 
field. 

In order to minimize the effect of spacecraft magnetic fields, the sensor 
portion of the magnetometer will be mounted at the end of a 3-6 n»ter extensible 
bo(mi to be provided by the spacecraft. Boom deflections and twist will be 
monitored by an attitude transfer system, provided by the spacecraft, which 
will taake available an orthogonal transformation between the sensor coordinate 
system and the coordinate system of a stellar attitude determination system 
within the spacecraft body. If required for magnetic cleanliness, the 
electronics for the magnetometer can also be mounted within the spacecraft 
body, and connected to the magnetcwieter by cabling routed along the boom 
structure. 

The magnetometer sensor shall have defined a geometric coordinate system, 
and means shall be provided through optically flat mirrors, or other optical 
devices, for determination of two orthogonal axes of the geof^etric coordinate 
system to one (1) arc second accuracy. The magnetometer shall measure projec- 
tions of the ambient field along three orthogonal magnetic axes, vjhose nominal 
directions within the geometric coordinate system shall be specified to one (1) 
arc second precision. 
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2. Enviroroiental Constraints - The environmental constraints on the sagneto- 
8»ter are presented in section 4 of this specification. All performance 
specifications of section 3 shall be met after exposure to shock, vibration 
and acceleration tests. All performance specifications shall be i«t under 
exposure to the orbital thermal and under ground sinxjiation of that environment. 
It is expected that both active and passive thermal design of the magnetaneter 
sensor will be necessary to satisfy the requirements of section 3. 

3. Performance Specification 

3.1 Stability of Magnetic Axes - The angular deviation of any magretic axis 
from its nominal direction shall not exceed 5 arc seconds after exposure to 
shock, acceleration, and vibration testing in accordance with sectitH) 4 of 
this specification. The angular deviation of any magnetic axis fr«n its 
nominal direction shall not exceed 5 arc seconds during exposure to the thermal - 
vacuum environment specified in section 4. 

3.2 Orthogonality of Magnetic Axes - The angle between any tmo magnetic 
axes shall be 90® plus or minus 0.1®. 

3.3 Range - The instrument shall be capable of measuring field coim>onents 

-■9 

along any axis from -64000 gansna to plus 64000 gamma. (1 gamma equal to 10 
Tesla). 

3.4 Resolution - The resolution along any axis shall be plus or minus 1 
gamma or less. 

3.5 Measurement Bandwidth - The instruir.ent shall be capable of measurement 
in a bandwidth of up to 25 Hz, and shall be compatible with bandwidth 
limiting to 1 Hz. 
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3.6 Nois e - Fluctuation noise associated with any axis shall be less than 
0.1 gamma zero-to-peak when measured for 10 seconds in a 25 Hz bandwidth. 

3.7 Zero-Offset - The zero offset change associated with any axis shall 
be less than plus or minus 0.4 gamma over the expected range of sensor and 
electronics temperature. Long term stability of the zero offset shall be 
less than 0.4 ganma/year. 

3.8 Absolute Accuracy - The total instnxnental error from all sources in 

the determination of the field component along any sensor geometric axis shall 
be less than 5 gamma. This error budget includes the effects of magnetic 
axis stability, zero offset, fluctuation noise, quantization noise, sensitivity 
changes, and all other purely Instrumental sources of error. 

3.9 Power Consumption - The total power consumption for the raagnet©neter, 
including power for normal control of the sensor, shall be less than 6 watts. 

3*10 Weight - Total weight of the system shall be less than 5 kilograms, 
inclusive of the boom cable interconnecting sensor and electronics. 

3.11 Volume - The sensor structure shall be coritair«d within a 30 cm. 
diameter sphere, while the electronics package shall be contained within a 
volume 20 cm X 15 cm X 15 cm. 


4. Environmental Requirements 
4.1 Vibration - 


4.1.1 Sinusoidal Vibr ation 
Axis 


Frequency 
Range (Hz) 


All (X,Y.Z) 


10-18 

18-36 

36-150 

150-2500 


Level 

(0 to Peak) 

5.0 g* 

0.3 in (O.A.) 
20. g 
5.g 


Sweep Rate 
(Oct/Min) 


Limited to O.o in double amplitude (O.A.) 
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4.1.2 Randexn Vibration 

Frequency 

PSD Level 

Acceleration 

Duration 

Axis 

Range(Hz) 

(g2/Hz) 

(q-rms) 

(minutes) 


20-200 

0.045 




200-400 

+3dB/0ct 

12.9 

2/axis 


400-2000 

0.090 




4.2 Acceleration - The performance specification of section 3 shall be 
satisfied after non-operating exposure of acceleration of 22. 5G in the thrust 
axis and 6.0G in the lateral axis. The acceleration shall be simultaneously 
applied for a duration of three minutes. 

4.3 Storage Temperature - Performance specifications shall not be con^romised 
by non-operating storage for 6 hours at -50“C and 6 hours at plus 85®C. 

4.4 Thermal Shock - The performance specifications shall not be compromised 
by exposure to 5 cycles of thermal shock, a cycle consisting of: 

Step 1 - 1 hour storage at -50"C 

Step 2 - transfer to plus 85“ C in not more than 5 minutes 

Step 3-1 hour storage at plus 85“C 

Step 4 - transfer to -50“C in not more than 5 minutes 

4.5 Electronics Operating Temperature - The performance specifications shall 
not be compromised by operation of the electronics at all temperatures in the 
range of -10“C to +50“C while the sensor is at room temperature conditions. 

4.6 Thermal Vacuum Operation - The temperature of the sensor is to be 
actively and passively controlled by the contractor to insure temperatures 
compatible with the performance specifications of section 3. In orbit the 
sensor will he exposed to complete insolation and complete eclipse on a varying 



duty cycle basis. The extremes of the cycling are expected to be (a) a 
continuous insolation for the indefinite period of time and (b) sixty (60) 
minutes of insolation to thirty (30) minutes of eclipse. 

Confirmation of performance under these conditions will be performed by 
the government during acceptance tests of the magnetometer. The sensor will 
be operated at pressure of 1/2 micron or less in a ncn-magnetic thermal vacuum 
system and magnetic test facility in which the environment specified herein 
shall be simulated by use of a liquid nitrogen shroud and solar simulator. 
During these tests, the electronics will be operated at atmospheric pressure 
and in the temperature range -10®C to +50°C. 
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tso logon rap 

E?:TER tso PASSWORD: aaBzasi !2 

RHP LOGO?.' IH PROGRESS AT UrflUiBS OK f^AY 26 , 7977 

WELCOME TO UCLA-CCK TSO 

READY: alloc fi(taykeep) da(»ykcep) 

READY: speakoz 

TSO SPEAKEASY 3 MU 2 :«I 5 PH MAY 26 , 1977 
; iibindex 

A 8 HEAS CHECK COHPAR EO EQUAT 3 MACSOL MAGCAL HACCALl 

MAGSIM HAINPO KMATRIX MODEL ? 7 UTHUT fi REWSECALriEWSIH SEMCAL 

'■EKCALl SOLEQU THAT TRAMAT 


''dit i“ 3 gsiR! 

H'GSIM IS HOT DEFI? 7 ED 
HA?!UAL ??ODE 
:__kspt{pagsim) 

:_edit magsira 
EDIT CC?-!NA?fD MODE 

• m 

• m 


1 

EDITING HAGSIM 

1.0 P. *>008 AM 

2.0 I?7PUT HUYX,?^UZX 

3.0 INPUT HUXZ,KUYZ 

4.0 MUXX=-1+SQRT(7-(NUYX*»2+MUZX*»2)) 

4.7 N'J2Z=-l+SQRT(1-(MUXZ»*2-i-:{UYZ**2)) 

4 . 2 DELHU=yEC( : MUXX , MUYX , HU2 X ) 

4 . 3 DELKU=VECC : NUXZ , HUYZ , KU2Z) 

4.4 XUV=VEC(:1 0 0) 

4.5 YUV=VEC(:0 1 0) 

4.5 ZUV=VEC(:0 0 1) 

4.7 MUrXUV+DELKU 

4.8 ?7U=2UV+DELHU 

4.9 PRIf{T(DELHU,MU,DELHU,t.'U) 

5.0 R=MAT(3,3:) 

5.0 INPUT J7E PLEASE 

7.0 ?RI?7T(f<E) 

8.0 K=VEC(72:) 

9.0 FOR 1=1, HE 

10.0 PRINT (’DATA FOR EXPERIMENT #’,I) 

11.0 INPUT R MATRIX PLEASE 

11.1 PRIHT(R) 

12.0 iIP=R»NU 

12.1 PRIKT(flP) 

13.0 M(I) = IfIMER(r.U,NP) 

13.1 PRIHT(M) 

*14.0 EiiDLOOP I 



OBKaNASPAiaa 
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A3 


keptCaageall ) 

: edit aagcall 
E&IT CCHHAKD MODE 
:nist 


pEOlTlHG HAGCAL1 
1.00 PROGRAM 
1.10 HARGIHS(132) 

1.20 IHPUT HE THE NUMBER OF EXPERIMEBTS 


1.3C 
1.31 
1.50 
1.60 
2.00 
3.00 
A. 00 
9.00 
10.00 
11.00 
12.00 
13.00 
15.00 
16.00 
17.00 
18.00 

19.00 

20.00 
21.00 
22.00 
22.50 

24.00 

24.50 

25.00 
25.10 
25.20 
26.00 
23.00 
29.00 
30.00 
31.00 
32.00 

33.00 

34.00 

34.50 

35.00 

35.50 

36.00 

37.00 

38.00 
38. 10 

39.00 

45.00 
46. CO 

47.00 
4B.no 

40.00 


AA=MAT(HE,4: ) 

HE=VEC(NE:) 

CA=HE 
A=VEC(4: ) 

PRINT {‘ARE OLD DATA AVAILABLE?*) 

INPUT AHS (0=F4 1=T) i 

IF(AKS.GE.l) GO TO LOCI 
FOR 1=1, HE 

HENCEFORTH R IS OBJECTC *R* ,1) 

PRIHTC INPUT DATA FOR EXPERIMENT f *.I) 

INPUT MECI) THE NORMALIZED MAGNETOMETER MEASUREMENT 
ME(I) 

IHPUT R THE TRANSFORMATION MATRIX 
R 

EHDLOOP I 

IHPUT RMSERR (THE RMS ERROR IN THE MORHALIZED MAGNETOMETER MEASUREMENT; 
VARYI=MAT(HE,1:) 

VARYI(,1)=RHSERR»»2 
RMSERR ;VARYI 
$ 

INPUT HUY, MUZ (THE Y AND Z COMPONENT OF HU) 

IHPUT HUX.HUY (THE X AKD Y COMPONENTS OF NU) 

HUX=SQRT(1-HUY»*2-MUZ»*2) 

KUZ=SQRT(1-HUX»»2-KUY«»2) 

MU=VEC(:MUX,HUY,HUZ) 

NU=yEC(:NUX,NUY,HUZ) 

MU; MU 

LOCI: CONTINUE 
$ 

FOR 1=1, HE 

HENCEFORTH R IS OBJECT( ’ R * , I) 

CA(I)=IMMER(HU,R»NU) 

A1=H»HU 

A2=TRANSP0SE(P)»MU 

A( 1 )=A1(2)-Ai:i )»MUY/MUX 

A(2) = A1C3)-A1(1)»MUZ/?!UX 

A(3)=A2(1)-A2(3)*HUX/HUZ 

A(4)=A2(2)-A2(3)*NUY/KUZ 

AA(I,)=A 

EHDLOOP I 

AA imSQfAL PAGB V 

RES=ME-CA Of POOR QUAUTl 

ME;CA;RES 

S 

SETHULLd .OE-10) 

HENCEFORTH T IS TRANSPOSE 
SEIGENAMALYSI3 OF AA»T(AA) 


• • » 1 A 



52.00 

53.00 
- 5^#00 

55.00 
56. CO 

57.00 

53.00 

59.00 

60.00 
61.00 

52.00 

63.00 

64.00 

65.00 

66.00 
6T.00 

63.00 

69.00 

70.00 
71.00 
72. CO 

73.00 

74.00 

75.00 

76.00 

77.00 

73.00 

79.00 
SO. 00 

81.00 
82.00 
83.00 

84.00 

85.00 

86.00 
87.00 

83.00 

69.00 

90.00 
91.00 

92.00 
93. GO 

94.00 
94.10 

95.00 

96.00 

97.00 

93.00 

99.00 
100.00 
101.00 
101 . 10 
102.00 
103.00 

104.00 
104.10 

105.00 
106.00 
106.50 
107.00 

107.50 
103. CO 

103.50 


AUEVTsHUEV A3-2 

0P*O 

INPUT LOMLIK (LOWLIM IS THE LOWEST ACCEPTABLE EIGEMVALUE) 

LOWLIH 

WHERECA'JEV.LT. (LOWLIM)) AUEVT=0 
lUsLOCS(AUEVT) 

IF(IU(1).EQ. 1) CO TO LOCAU 
I=LOCS(.KOT.AUEVT) 

UP=U(,I0) 

UO=U(.I) 

LOCAU: CONTINUE 
$EICEKAHALYSIS OF T(AA)»AA 
AV=T(AA)»AA 
AVEV=EIGEHVALS(AV,V) 

AVEV 

AVEVTsAVEY 


VP=V 

WHEBE( AVEV. LT. LOWLIM) AVEVTsO 
IV=LOCS(AV£YT) 

IF(IV(1).EQ. 1) GO TO LOCAV 
IrLOCSC.KOT.AVEVT) 

VP=V(,IV) 

VO=V(,I) 

LOCAV: CONTINUE 

$SEE IF UP AND VP GIVE A POS DEF DIAG HAT 

L=T(UP)»AA»VP 

DE=DIACELS(L) 

WHERE (DE.GT.O) DE=0 
IFXrLOCS(DE) 

$CHANC£ SIGN OF COLUMNS OF UP CCESESPOHDIHG TO IIEGATIVE EIGENVALUES 
HDO=;iOELS(IFX) 

IFCIFX(I).EQ.O) GO TO BELOW 
FOR 1=1, KDO 

UPC ,IFX(I))=-UP{ ,IFX(D) 

EMDLOOP I 
BELOW: CONTI HUE 
L=T(UP)*AA*VP 
DIACELS(L) 

IMVAA=VP*(1/L)ST(UP) 

X=IHVAA»RES 

NR=MOROWS(INVAA) 

HC=NOCOLS(IflVAA) 

H2=HAT(NR,NC: ) 

FOR 1=1, NR 
FOR J=1,NC 

H2(I,J)=IHVAA(I, J)*»2 
ENDLOOP J 
ENDLOOP I 
yARXK=H2»VARYI 


SQRT(VARXK) 

IF(IU(1).EQ.l) GO TO HERE 

VAR=(T(U0)»RES)»»2 

LRESSQ=RES»»2 

VAR;LRESSQ;VAR/LRESSQ 

HERErCONTIHUE 

X 


OF POOR QUAWinr 


MUY=X(1)+MUY 

MUZ=X{2)+MUZ 

?JUX=X(3)+NUX 

f!UY=X(4)+?!UY 

rtUX=3QRT(1-MUY*»2-MUZ**2) 
NUZ=SQBT(1-HUX**2-?IUY»*2) 
I-IU^EC ( rMUX.MUY.MUZ) 


1 
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113.00 .MU;HU»»2 

m.OO P8I»T(»ARE THESE VALUES GOOD EMOUGH?*) 
-Il 5 ,e 0 IHPOT Ar?S (0=F & UT) 

116.00 IF(AKS) GO TO PROEHD 

117.00 GO TO LOCI 

118.00 PROEHD rCOIITIHUE 

119.00 AA;IHYAA;U;V;TKVAA*AA;AA»IKVAA;H2 
119.10 IFCIO(l).HE.l) PRIHT(V0.UQ,T(U0)*RES) 

•120.00 PBISTCMOSaAL EHD OF PROGRAM HACCAL*) 
1 % 

end 


P 


Mter apl or tso 
tso 

^0/91 temporarily not available 

6 



J^endix M 

• 

end 

MANUAL NODE 
:_edlt newsla 
KEWSIM IS HOT DEFINED 
MAlfUAL NODE 
;_k«pt(newsitn) 

: edit newsim 
El?IT COMMAND MODE 
:S 


1 

EDITING NEWSIH 

1.0 PSOGSAM 

2.0 INPUT MUT A HAT WITH ROWS BEING UNIT VEC OF SEN IN SEN CCOR 

3.0 INPUT HUT A MAT WITH ROWS BEING UNIT VEC OF COIL IN GEO COOR 
A.O HENCEFORTH T IS TRAKSP 

5-0 MU=T(KUT);HU=T(NUT) 

6.0 HU;HU 

7.0 SOBTAIM LIST DEFIMIrlG NOH OF SEN ORIENT AND TRAKS MAT 

8.0 KEPTLIST(THAT) 

9.0 R0HS=9*SE 

10.0 B=VEC(ROWS:) 

11.0 $ 

12.0 $LOOP TO CREATE MEASUREMENT HAT FOR CAL MAT BC=I 

13.0 L=1 

U.O FOR Kr1,KE 

14.1 PRINT (* EXPERIMENT # =', K) 

15.0 HENCEFORTH R IS OBJECTC ’ R ’ ,K) 

16. C BM=MUT»R»HU;BM 

17.0 5 

15.0 $ PUT COLUMNS OF BM INTO COLUMN VEC OF HEASUREMEHTS 

19.0 FOR J=l,3 

20.0 B(L)=BM(,J) 

21.0 L=L+3 

22.0 ENDLOOP J;EHDLOOP K 

23.0 SPRINT AND SAVE THIS VECTOR 

24.0 TMAT=NAMELIST(NE,R1,R2 R3 R4 R5 R6,R? ,B) 

25.0 KEEPLIST(TMAT) 

»26.0 B 

• CP 

• 


end 

MANUAL flODE 


OIOGINAL 
OF POOR OUALITY 
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« ^p^r.d1x AS 

EOITIiJG MAIHPO 

1 PROGRAM 

2 $ FROM HERE TO ITLOOP: IS PROGRAM IMITIALI2ATIOH 

3 DELHU=MAT(3.3:);DELNU=DELHU 

tt PRIMK* INITIALIZE PROGRAM TO CALCULATE SENSOR AND COIL COSINE 

5 INPUT HUT, HUT (ROWS ARE GUESSES OF SENSOR AND COIL UNIT VECTC 

6 HENCEFORTH T IS TRANSPOSE 

7 MU=T(MUT);fiU=T(NUT) 

S PRIHTCTHE INITIAL MODEL ON INPUT WAS') 

9 MU;NU 
10 $ 

11 $ OBTAIN MEASURED TRANSFORMATION MATRICIES 

12 RRrHAT(21,3:) 

13 KFPT(TRAMAT) 

1A LCADOATA(RR,TRAHAT) 

15 PRIHTCTHE INPUT TRANSFORxMATIOH HATRICIES ARE') 

16 RR 

17 $ 

18 $ 03ATAIM MEASURED FIELD VALUES 

19 KE?TLIST(BMEAS) 

20 B 

21 $ 

22 $ FREE AS MUCH SPACE AS POSSIBLE 

23 FREE (MUT, NUT, TRAHAT) 

24 $ 

25 $ FIND OUT HOW MANY SENSOR ORIENTATIONS MEASURED 

26 INPUT NE (THE NUMBER OF SENSOR ORIENTATIONS) 

27 R0WS=9*NE 

28 S 

29 S COME HERE FOR ITERATION OF SOLUTION 

30 LST0P=3 

31 LM?=1 

32 ITLO : 

33 $ OBTAIN AH EXECUTE PROGRAM TO CREATE EQUATIONS FOR THIS MODEL 

34 PRIiTTCBEGINMISG ITERATION NUMBER ',LMP) 

35 KEPT(EQUATS) 

36 EXECUTE (ECU ATS) 

37 FREE(EQUATS) 

33 $ OBTAIN AND EXECUTE PROGRAM TOS SOLVE EQUATIONS 

39 $:<EPT(SOLEQU) 

40 $EXECUTE(SOLEQU) 

41 $F8E£(SCL£QU) 

42 KEPT MACSOL 

43 EXECUTE MACSOL 

44 FREE MACSOL 

45 LMP=LMP+1 

46 IF(LMP.LE.LSTOP) GO TO ITLOOF 

47 MU =M Ui-DE LHU ; HU =N?J+DELMU 
4*3 DE LflU * DELNU ■ MU * NU 

49 PRliiTCIfJITIAL’HODEL HAS BEEN IMPROVED LSTOP TIMES ' ,LSTOP) 


(/I tc 



edit equats 
EDIT COMMAND MODE 

. C I * 


1 1st 


EDITING EQUATS 

1 PROG3AM EQUATS 

2 $ UPDATE MU AMD HU WITH PREVIOUSLY CALCULATED CORRECTIOMS 


3 MU=HU*DELMU;NU=?!U+DELMU 
A DELMU;DELKU 

5 free(delmu;delnu) 

6 $ CALCULATE DIAGONAL ELEMENTS OF HU AND HU USING OFF DIAG ELE 

7 HU(1 , 1)=SQRT(1-KU(2,1)*»2-HU(3, 1)**2) 

B ;;U(1 , 1)=SQRT(1-NU(2,1)»»2-NU(3,1)**2) 

9 KU(2,2)=SQRT(1-MU(1,2)»»2-MU(3,2)**2) 

10 MU(2,2)=SQRT<1-HU(1,2)»»2-NU{3,2)**2) 

11 MU(3,3)=SQRT(1-HU(1,3)**2-HU(2,3)**2) 

12 KU(3,3)=SQRT(1-NU(1,3)**2-HU(2,3)*»2) 

13 PRINT! »MU AHD MU AS CORRECTED TO GIVE UHIT COL NORMALS ARE') 

14 ML«;NU 

15 SCREATE UHIT VECTOR COHSTRAIKT MATRIX 

16 V1=NU(,1)/KU(1,1) 

17 V2=HU(,2)/NU(2,2) 

18 V3=HU(,3)/HU(3,3) 

19 V4=MU(,1)/MU(1,1) 

• 20 V5=MU(,2)/HU(2,2) 

21 V6=MU(,3)/MU(3,3) 

22 $D0 LOOP FOR EACH ORIENTATION 

23 PRINT! 'BEGINNING LOOP FOR EACH SET OF SENSOR ORIENTATIONS') 

24 J= 1 5 9 10 14 18 

25 L = 1 

26 A =MA7 ! ROWS , 1 2 : ) ; BR =VEC ! ROWS : ) 

27 FOR :<=1,NE 

28 AA=MAT!9,18: );NN=MAT!9,18: );RNX=MAT!7,3: );BM0D=VEC!9: ) 

29 RNX=HAT!7,3:) 


30 $ PUT PROPER TRANSFORMATION MATRIX III R 

31 IT=3*!X-1)+1 

32 I3A=IT,IT+1,IT+2 

33 R=RR!ISA,) 

34 SCALCULATE MODEL MATRICIES 

35 C1=T!MU)»R;C2=R*HU;BMM=T!MU)*R»NU 

36 RfiXCl )=C2! ,1) 

37 RNX(4)=C2! ,2) 

38 PNX(7)=C2! ,3) 

39 S SET UP MODEL MEASUREMENT VECTOR 

40 3M0D(1 )=BHM( , 1) 

41 BMOD(4)=BHM( ,2) 

42 BM0D(7)=BHK! ,3) 

43 I3A = IfITEGERS!L,L+8) 

44 B3=B!ISA) 

45 PR(L)=BB-BMOD 

46 33ET UP LEFT HALF OF COEFFECIENT MATRIX 




T 


47 AA(1,1)=C1 

48 A.A(4,4)=C1 

49 AA(7,7)=C1 

50 $35T UP RIGHT HALF OF COEFFECIENT HAATRIX 

51 AA( 1 , 10) = R?iX 

52 AA(2,13)=RHX 

53 A\( 3, i5) = Rrix 

54 i:?i( 1 , 1) = CCLARRAY(9:AA( , 1 ) ) *ROWARRAY( 3 : VI ) 

55 ,4)=C0LARRAY(9:AA( , 5 ) ) * ROWA.RRA Y( 3 : V2 ) 

55 ::"0 ,7)=COLAPP.AY(9:AA( ,9))"RCWARRAY(3:V3) 

57 r;:i ( • , 10) = CCLA RRAY(9:A A( , 10 ) )^RQWARRAY(3 : V4 ) 



rtNl I , l3i=UULAKKAUy :AAl , lOii 
60 AAsAA-HFAMCHN) 

6! /A=EUMC0LSCftA, J) 
62^A(L,1)sAA 
63 L=L+9 
6A E.'iDLOOP K 


KUWAHKAi Ij: tfb) 



65 BR 

65 FREECV1 ,V2,V3,Vi|,V5,V6,J,K,AA,MK,IT,I3A,R,C1,C2,BMM,RNX,BMCD,BB) 

»67 PaiMK ’PROGRAM EQUATS HAS COMPLETED CALCULATION OF EQUATIOHS FOR THIS 
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end 

.MAHUt'L NODE 
•:_iccV'- ^ solequ) 

: edit solequ 
EDIT COMMAND MODE 
: SI ist 


EDITING SOLEQU . 

1 PEOGRAK 

2 $ INITIALIZE ARRAYS FOR SOLVING SET OF EQUATIONS 

3 X = ’. t;C(ROWS: ) 

4 T=1.0E-6 

5 S=ARKAY(l2:) 

6 X=SIMEQUAT(A,BR,T;S) 

7 0 

8 $ SET UP CORRECTION MATRICIES 

9 OELNU=TRANSP(MAT(3,3:0 X(l) X(2) X(3) 0 X(4) X(5) X(6) 0 )) 

10 0ELMU=TRAI!SP(MAT(3,3:0 X(7) X(8) X(9) 0 X(10) X(11) X(12) 0 )) 

11 FREE(X,T,S,A,BR) 

12 $ 

»13 PRINT! ’PROGRAM SOLEQU HAS COMPLETED SOLUTION OF EQUATIONS FOR THIS MODEL’) 



